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Strahlungsenergie und Erdabsorption 
bei Dipolantennen 


Von A. Sommerfeld und F. Renner 
(Mit 5 Abbildungen) 


$1. Einleitung und Übersicht 

Die Strahlungsquelle sei ein monochromatisch schwingender 
Hertzscher Dipol von vertikaler oder horizontaler Achsenrichtung. 
Die Erdoberfläche werde als eben angenommen. Die Beschaffenheit 
der Erde werde durch die zunächst beliebigen, aber innerhalb der 
Erde festen Konstanten « und o.gekennzeichnet. (u sei gleich 1 
gesetzt) Indem wir wie Hertz in Gaussschen (nicht-rationellen) 
Einheiten rechnen, nennen wir die Wellenzahlen im Metlium I (Luft) 
und II (Erde) 

(1) 

Der Abstand des Dipols von der Erdoberfläche heiße h. Wir 
werden die Abhängigkeit der Strahlungsenergie usw. von h und 
ihren von der Erde absorbierten Anteil durch verhältnismäßig ein- 
fache Formeln festlegen und durch Abbildungen yeranschaulichen. 

Die allgemeine und strenge Darstellung des Feldes ist bekannt. 
Im Falle des: horizontalen Dipols geht sie auf Hörschelmann') 
zurück. Besonders das Feld des vertikalen Dipols ist vielfach dis- 
kutiert, in der Absicht, die in der strengen Darstellung vorkommen- 
den Besselschen Funktionen zu eliminieren oder durch elementare 
Funktionen zu approximieren. Unserer Aufgabe sind aber gerade 
die strengen Formeln angepaßt. 

Es handelt sich bei uns um die Integration des Energieflusses 
durch eine den Dipol umgebende Fläehe bzw. durch ein Flächen- 
paar. Bei vollkommen leitender Erde liegt es nahe, als solches 
die unendlich ferne Kugel zusammen mit der spiegelnden Erd- 
oberfläche zu benutzen (Abraham 1899). Bei beliebiger Erde wird 
dies aber wegen Komplikation der Formeln untunlich. Vielmehr 
ist es angezeigt und der Symmetrie des Problems angepaßt, bei 
Berechnung der Gesamtausstrahlung ein den Dipol umschließendes 

1) H. von Hérschelmann, Münchener Diss. 1911; Jahrb. d. drahtl. 
Telegraphie u. Telephonie 5. S. 14. 1911. 
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Ebenenpaar zu benutzen, vgl. Abb. 1. Zur Berechnung der Erd- 
absorption dient dann die untere der beiden Ebenen. Da der Energie- 
fluß durch das Produkt zweier Feldkomponenten gegeben ist, deren 


4 


Abb. 1. Der Ursprung 0 des rechtwinkligen xyz-Systems liege senkrecht 
unter der Antenne A. Beim horizontalen Dipol weise x in die Achsenrichtung 


desselben. Daneben benutzen wir Zylinderkoordinaten r, 9,2. Der Winkel 9. 


werde von der x-Achse aus gezählt. R sei der Abstand des Aufpunktes P vom 

Ort der Antenne: R = Vr' + (¢ — h)*. In der Abbildung ist auch das Spiegel- 

bild A’ des Antennenortes A eingezeichnet. Der Abstand des Aufpunktes von 
diesem ist R’ = Yr* + (z + h)* 


jede mit einer Besselfunktion J, (Ar) bzw. J, (lr) behaftet ist (4 und 1 
sind weiter unten zu erklärende Integrationsvariablen), so kommt es 
durchweg auf Integrale der Form an 


(2) I ft 


In der Tat en wir zeigen, daß die beim "horizontalen Dipol 
ursprünglich auftretenden, etwas komplizierteren Integrale sich auf 
die Form (2) zurückführen lassen. Beim vertikalen Dipol tritt direkt 
das Integral (2) mit m = 1 auf. Es gilt nun einfach: 


--b=A, 


so zwar, daB bei endlichem. oder beliebig kleinem 7 gilt: 


7 
(3a) „ma=ı. 
I hat also den Charakter einer „ö-Funktion“. 

Auf Grund von (3) und (3a) läßt sich nun die Integration des 
Energieflusses über die Ebenen z=h+s ausführen. Dabei ver- 
schwinden die Besselschen’ Funktionen aus unserer Formeln und 
reduzieren sich die ursprünglich dreifachen Integrale nach r, ! und A 
auf einfache Integrale. 

Wir nennen I das Orthogonalitätsintegral der Besselschen 
Funktionen. Dieses gestattet, wie immer, die Darstellung von will- 
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kürlichen Funktionen»), Übrigens geht die Bedingung I = 0 aus 
der Orthogonalitätsbedingung der Kugelfunktionen 


+1 
-1 


durch einen bekannten Grenzübergang hervor, der bei unserem 
Problem genau dem Übergang von der Erdkugel zu der ebenen Erde 
entspricht. Es war dieser Grenzübergang, der uns die Integration 
über die Ebenen der Abb. 1 nahelegte. 

Dieselbe Integrationsmethode ist schon früher von K.F.Niessen?) 
durchgeführt worden. Auf einige Einwände?) gegen seine Ergebnisse 
kommen wir in §§ 4 und 5 zurück. Wir glauben, daß unsere Be- 
handlung einfacher und vollständiger ist als die Niessensche. 

In 8 3 wenden wir unsere Integrationsmethode auf die gesamte 
Ausstrahlung S des vertikalen bzw. horizontalen Dipols an. Bedeutet — 
S, den gesamten Energiefluß durch die Ebene z=h+e in der 
positiven z-Richtung, S_ denjenigen durch die Ebene z=h-, 
ebenfalls in der positiven z-Richtung, so ist S gegeben durch 


„Zuvor müssen wir in § 2 die allgemeine Felddarstellung für 
beide Fälle bereitstellen, weil diese in der Literatur vielfach‘) nur 
im Grenzfall k = 0 angegeben ist, in dem unsere Integrations- 
methode hinfällig werden würde. Beim vertikalen Dipol wird das 
Feld abgeleitet aus einer Hertzschen Funktion II = II,, beim 
horizontalen Dipol (Achsenrichtung z) muß man bekanntlich einen 


Hertzschen Vektor i] = (I7,, 0, IT,) benützen. Das Gegenstück im 
Fall von Kugelproblemen allgemeinen Charakters ist das Auftreten 
zweier Teilfelder und zweier zugehöriger Potentiale’). 

In § 4 wird die Berechnung der Ausstrahlung S vervollständigt. 
Es wird nämlich neben dem Gliede nullter Ordnung in k,/|k,|, 
welches schon in $ 3 abgespalten war und die Ausstrahlung bei 
unendlich gut leitender Erde darstellt, das Korrektionsglied erster 


1) Vgl. z.B. Frank-Mises, 2. Aufl. 1935, 8. 921, im folgenden mit 
Fr.-M. zitiert. 

2) K. F. Niessen, Ann.d. Phys. [5] 22. 8. 162. 1935; 24. S. 31. 1935; ° 
32. S. 444. 1938. 

3) In einer Korrespondenz mit Herrn Niessen hat sich eine völlige 
Übereinstimmung der beiderseitigen Standpunkte ergeben, wofür wir ihm auch 
an dieser Stelle danken möchten. 

4) Z.B. Fr.-M. Kap. XXIII. §3 1 u. 2. 

5) Fr.-M. Kap. XX. S. 873. 
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Ordnung ermittelt, welches z. B. für Seewasser auch bei kürzeren 
Wellenlängen, z.B. 4 = 40 m, den Sachverhalt hinreichend genau 
darstellt. In Zusatz 6 wird auch die nächste Näherung kurz 
diskutiert. 

85 behandelt die Frage der Erdabsorption, die wir mit Says, 
bezeichnen. Man berechnet diese als Ausstrahlung durch die 
untere Ebene z= h—e des in Abb. 1 dargestellten Ebenenpaares, 
und zwar in Richtung der. negativen z-Achse. Es wird also nach 
der vorangehenden Verabredung über die Bedeutung von Sx: 

San. = — S_. 
Da der unendlich gut leitende Erdboden ein vollkommener Reflektor 
ist, verschwindet S,y.. im Limes k, —» oo. Das Hauptglied von S ays, 
wird daher von der Ordnung k,/|k,|. Ergänzungen ‚zur Berechnung 

dieses und der höheren Glieder bringt Zusatz 6. 

In § 6 werden die beim vertikalen und horizontalen Dipol 
gewonnenen Resultate diskutiert, verglichen und graphisch erläutert. 

Der Schlußparagraph betrifft die Frage: Wie ist die Stärke 
unseres Dipols zu wählen, um diesen einer Antenne von gegebener 
Länge und Stromstärke bei .bestimmter Stromverteilung anzugleichen. 
Hier wird auch eine neue Formel für den sogenannten Strahlungs- 
widerstand entwickelt, welche die früher hierfür gegebenen Formeln 
mit Rücksicht auf die Bodenbeschaffenheit und die Lage der Antenne 
erweitert. 


$2. Allgemeine Formeln zur Darstellung des Feldes einer Dıpolantenne 


Bei Unterdrückung des Zeitfaktors, den wirin der Form exp. (—imt) 
annehmen, gilt für das Medium I (Luft) 


> -> 
AU +k?U=0, 
(4) > > 
ik rot. 
Beim vertikalen Dipol ist 
&=(E,E), 9=H 
also wegen (4), wenn wir hier und im folgenden nur diejenigen 
Komponenten anschreiben, die wir direkt brauchen werden’): 


0° II, OT, 
©) He = 
1) Die Formeln für IZ, € und § im Medium II (sowie für E, im Medium I) 
sind zwar unentbehrliche Bestandteile der vollständigen Lösung, kommen aber 
für die folgenden Rechnungen nicht direkt in Betracht. Sie werden in Zusatz 1 
“am Ende der Arbeit) kurz zusammengestellt. 
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Um /7, zu finden, gehe man von der primären Erregung aus 
(6) Trim. = A 


R 
und bemerke, daß 


Partikulärlösungen von (4) in I bzw. II sind, die der Symmetrie 
unseres Problems und der Forderung des Verschwindens für z = +00 
entsprechen. Indem man ‘mit je einer frei wählbaren Funktion f(2) 
multipliziert und über A von 0 bis oo integriert, erhält man einen hin- 
reichend allgemeinen Ansatz für die „sekundäre Erregung“, um mit 
diesem die Grenzbedingungen für z = 0 (Stetigkeit von E, und H, beim 
Übergang zum Medium II) erfüllen zu können. Man findet für das 


Medium I: 

mit 

(8a) M = k,*u, — k,? 


J ohne Index bedeutet hier ynd im folgenden die nullte Besselsche 
Funktion. In (8) gilt das obere Vorzeichen für z>h, das uutere 
für 0<2z<h. Man beachte, daß für z= h, wie es sein muß, die 
beiden in (8) zusammengefaßten Darstellungen ineinander übergehen. 
In (8) ist das erste Glied rechts identisch mit der primären 
Erregung (6). Das zweite Glied geht für k,—-> oo über in die 
spiegelbildliche Erregung (vgl. Abb. 1) 
R’ 
Nach (5) erhalten wir aus (8) Integraldarstellungen für E und H,. 
Wir schreiben dieselben —' für die beiden Ebenen z=h+e an: 


fra rf, (Marar, 
(9) 
(ar) f,(a)arda. 
0 
Hier bedeuten: 


(9a) 
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In den zweiten Gliedern rechts von (9a) heißt der Exponent ur- 
sprünglich — u, (2h + e). Man kann aber s neben 2h streichen und 
braucht es nur da beizubehalten, wo es allein auftritt, d.h. in den ersten 
Gliedern rechts. 

“Beim horizontalen Dipol ist (vgl, die Einleitung): 


if =(i,, M1). 
Aus (4) ergibt sich 


0° II, 
(10) E, = cong + h,? cosy IE, 


sing 9I, 


= 
E, k,*sin IT, Dr 


(10a) 


Die primäre Erregung, die nur in /7, vorkommt, setze man an als 
(11) 


Neben (7). braucht man hier als Partikularlösungen, die der Sym- 
metrie des horizontalen Dipols angepaßt sind: 
cos J, (Arjetm™=, 


cos J,(Ar)etm:, 


Statt J, kann man bekanntlich auch — J,’ schreiben. Aus solchen 
Partikularlösungen baut man Integrale auf, die man durch Wahl der 
in sie eingehenden verfügbaren Funktionen von A den Grenz- 
bedingungen für z = 0 auch hier anpassen kann. Man findet so für 
das Medium I: 


(13) 


Wegen der Bedeutung des doppelten Vorzeichens ‘sustain man die 
obige an Gl. (8) anschließende Bemerkung, ebenso wegen des stetigen 
Übergangs bei z=h. N ist dieselbe Abkürzung wie in (8a). 

Aus den Gl.(10) und (10a) folgt dann für die beiden Ebenen 


= 
—H, = c08 9 — =", 
? 
i 1 01, + sir ol, 
4 
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¥. 


Fas = [ann @edr— fans, rar, 
+ (a) a? da + Mada, 
(15) 
0 


Die Bedeutung der hier eingeführten Abkürzungen f,, f,,... ist, 
wenn man [vgl. oben bei(9a)] sogleich e gegen 2h streicht, die folgende 


etm 


2 1 
(16) 


etn), 


(17) ‘ N 


'„$3. Berechnung des totalen Energieflusses. 
Reduktion der dreifachen auf einfache Integrale 


Da der Energiefluß in Gaussschen Einheiten durch [9] 


gegeben ist, kommt es in unserem Falle (vgl. Abb. 1) auf die Be- 
rechnung von 

(18) [EG] = 6.5, 8,6, 

an. Hier ist z. B. unter & der reelle Teil der (vgl. den Anfang 
von § 2) mit dem Zeitfaktor exp.(— iwt) versehenen komplexen 
Größe E, gemeint, also der Ausdruck 


entsprechend ist 
Eo, = (Be e* (He + Hy 


Der zeitliche Mittelwert, auf den es uns bei der Berechnung der 
Ausstrahlung allein ankommt, beträgt: 


(19) {E} H,}. 


8,5” 
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Re bedeutet hier und im folgenden „Reeller Teil von“. Ebenso gilt 
natürlich: 


er 1 
(19b) E,0, Re 
Wir berechnen daraufhin für den vertikalen Dipol den gesamten 


zeitlich gemittelten Energiestrom S; durch die Ebenen z=h+e, 
vgl. Abb. 1. Nach (18) ist wegen €, =§,=0: 


8. 9, (Re [EH, rdr 


oder mit den Ausdrücken (9) für E, und H, he mit fr ty 


Hier muBten wir die Integrationsvariable i in der zweiten Gl. (9) um- 


“benennen (l statt 2), um sie vonder in un ersten Gl. (9) zu unter- 
scheiden. 

Jetzt verwerten wir die Orthogonalitätsbedingungen (3), (3a). Sie 
reduzieren das vorstehende dreifache Integral auf 


(208) f FOL @: 
0 


Beim Einsetzen aus (9a) entstehen in dem Produkt f,*(a) f, (2) zu- 


nächst vier Glieder. Bilden: wir aber den Energiefluß durch das 


Ebenenpaar z =h +, nämlich 
(21) S=8,-S_, 


so°bleiben nur zwei Glieder übrig, und man erhält aus (20): 


Hier Babes wir in dem zweiten Glied, wie oben in (9a), & gegen k 
vernachlässigt. 

Das erste Glied läßt. sich sofort auswerten, wenh ‘man die 
Realitätsverhältnisse berücksichtigt. Es ist.nämlich nach (7) u, reell 
für 4>k,. Da natürlich auch u,* + u, reell ist, so kommt nach 


Multiplikation mit — i und Bildung des reellen Teils die Integration 
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von k, bis 20 in Fortfall. Die verbleibende Integration von 0 bis k, 
läßt sich, zumal da hier #,* + a, = 0 ist, sofort ausführen und gibt?) 


0 


Von dem zweiten Gliede des Integrals in (22) wollen wir den- 
jenigen Bestandteil vorwegnehmen, der beim Guenzübergang k, —> co 
nicht verschwindet. Wegen der Werte (8a) fir M und "N wird 


=—1fürk, —o. Wir ersetzen daher in (22) durch 


1- (+ + 1); dadurch dutsiheh ‘die beiden Integrale 


(22b) J Js ane 


und 
a 
Das in Rede stehende zweite Glied in (22) wird dann gleich 
Re(i J) + Rö(iK) 
Bezüglich der Realitätsverhältnisse gilt fir J dasselbe wie beim 
ersten Gliede. Die Integration darf also auch hier auf den Be- 


reich 0<A<k, beschränkt werden und läßt sich dann leicht aus- 
führen. Macht man die Substitution 


(224) p, =- 


1) In den grundlegenden Formeln (8) und:(13) ist die Zweideutigkeit 
von u, und u, allgemein so festgelegt, daß der Realteil bei beiden immer 
positiv sein soll. Um diese Vorschrift für 0 <A <k, zu ergänzen, .wo u, rein 
imaginär ist, geht.man so vor: Man nimmt zunächst schwache Absorption in 
Luft an. &,* liegt dann ebenso wie &,* im I. Quadranten, A? — k,®. im 
III. Qundrinten, u, also im IV. oder II. Weil aber der Realteil positiv sein 
soll, muß man den im IV. Quadranten liegenden Wert nehmen; daraus folgt, 
daß man auch bei fehlender Absorption das alsdann rein imaginäre u, mit 
negativem Vorzeichen rechnen muß. . 

Wir bemerken noch, daß entsprechend auch u, und im besonderen Y — k,?, 
(das ist der Wert von mw, für 2= 0) im IV. Quadranten liegen soll. Daraus 
folgt, daß für den in (34) einzuführenden Winkel 4 gilt 
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so transformiert sich J (mit k, als oberer Grenze) in 
1 
(22e) | im ik? fete (u — 1) du. 
0 
Hier haben wir die Abkürzung benutzt 
(22 f) C= 2k,h. 
(22e) gibt ausgerechnet 
(22g) J=k, (1+ i=) : 
Mit Rücksicht auf (22a), (22g) und die Bedeutung von K in (22c) 
geht (22) über in 
28 2 2 /sind. 

(23) — cost) + Re K}. 

Andererseits haben wir beim horizontalen Dipol wegen (19a,b) 


[EG], = + Re — 


Zu Sz übergehend führen wir zuerst die Integration nach 9 aus, 
die wegen (14) und (15) einen Faktor a liefert, also 


S,= Refrdr H,— E,H,} 
0 


und 


(24) = Re - 


Die hier angedeutete Klammer besteht nach Ausmultiplizieren der 
Ausdrücke (14) und (15) aus acht Gliedern, welche außer Bessel- 
schen Bestandteilen die Produkte 


enthalten. Diese sind aber nach (16) und (17) in S, und S_ die 
gleichen, mit Ausnahme von denjenigen, die den Faktor f, enthalten, 
wobei auch von f, nur das erste mit + behaftete Glied in Betracht 
kommt, das wir in der hier zu benutzenden Integrationsvariablen | 
schreiben werden: 

em, om, = — 
Gehen wir also, wie in (21) zu S=S, — S_ über, so vereinfacht 
sich die eckige Klammer in (24) erheblich, nämlich auf 


2e- «m J (Ir) [a J’ @+ 
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Hier lassen sich noch die beiden ersten Terme der { } nach der 
Besselschen Differentialgleichung zusammenziehen zu 

— (An f,*(Q). 
Somit folgt aus (24) 


ER Refi rar fran 


(25) 


Jetzt verwerten wir wieder die Orthogonalitätsgieichungen (3), (3a) 
und erhalten statt des vorstehenden dreifachen Integrals einfach: 


0 


Indem wir dies in (25) eintragen und hinter dem Zeichen Re kon- 
sequent — i durch + 7% ersetzen, was bequem ist, haben wir 


(26) - - if (A) + 2f, 
0 


Nach (16) ist 
(4) + 2f, (4) = 


— 
My, 


(26a) | 
1 - 21° 
Hier wollen wir von dem zweiten Teil denjenigen Bestandteil ab- 
spalten, der beim Grenzübergang k, —> 00 nicht verschwindet, das 


ist — — (2k,?— A*)exp.(— 2hu,). Wir erhalten so 


(27) (@) + 2f, = e-2hm) + 2F (A) e-2hm, 
Dabei bedeutet F(A) zunächst die folgende Abkürzung 


2k? — —u AP — 


Wenn man hier auf der rechten Seite den gemeinsamen Nenner 
(4, + %)N herstellt und dann Zähler und Nenner mit yp, —p, 
multipliziert, so ergibt sich einfacher: 


(27b) =k? 


V 
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Einsetzen von (27) in (26) liefert nun, wenn wir wieder e neben 2h 


vernachlässigen: 


(28) 
+ 24 


Das erste Integral läßt sich unter Berücksichtigung der Realitäts- 
verhältnisse (die obere Grenze ist wieder wie oben durch? k, zu 
ersetzen) leicht auswerten. Mit der Substitution (22d) und der. Ab- 
kiirzung (22f) witd es 


i + u2)(1 — = |} (1 
(28) ap daraufhin über ir 
mit der Abkiirzung 
(30) Le - 
0 


In (23) und (29) haben wir die vollständige Darstellung der Aus- 
strahlung eines vertikalen und horizontalen Dipols vor uns. Die 
ausgerechneten ersten Glieder hätten, da sie sich auf eine unendlich 
gut leitende Erde beziehen, auch ohne unseren analytischen Apparat 
gefunden werden können. Wesentlich sind für uns die zunächst 
noch in Integralform geschriebenen letzten Glieder, da sie ‘den 
Einfluß der endlichen Leitfähigkeit beschreiben. 


§ 4. Untersuchung der Korrektionsglieder 


Wir nennen die jetzt zu. betrachtenden Glieder Korrektions- 
glieder, weil sie von der Ordnung k ,/k, oder von höherer Ordnung 
klein sind. Zur vorläufigen Orientierung sei bemerkt, daß für die 
ne von Seewasser und bei einer Wellenlänge von 
etwa 40 m = ~ 5 ist. 

Wir wollen uns im Text auf die Glieder erster Ordnung in k,/k, 
beschränken. Wegen der Glieder zweiter Ordnung vgl. Zusatz 6. 
Wir haben zu untersuchen beim: vertikalen Dipol das Integral K, 
vgl. (22c) und (23), beim horizontalen uae das Integral L, vgl. (27b) 
und (30). 


. 
/ 
: 
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Da wir zunächst das Korrektionsglied nur in erster Annähe- 
rung bekommen wollen, können wir beim horizontalen Dipol die 
folgenden Vernachlässigungen machen. Im Zähler von (27b) kann 
man A? gegen 2u, u, streichen, obwohl A von 0 bis o geht, weil 
nämlich große Werte von 4 wegen des Faktors exp. (— 2hy,) bei 
der Integration unterdrückt werden. Im Nenner wollen wir 


(31) gegen k,*p, 
streichen. Die beiden Vernachlässigungen sind sicher falsch an der 
Stelle A=k,, also für u, =0. Daß aber die Umgebung dieser 
Stelle zum Integral keinen wesentlichen Beitrag liefert, zeigen die 
Rechnungen in Zusatz 6a. 

Nach diesen Verabredungen are wir 


Wir ersetzen noch konsequenterweise 
(32) #, durch Y—k,?. 
Damit erhalten wir aus (30) 
2k? 
Das entsprechende allgemeine Integral gibt ausgerechnet: 


(826) — + 


Wir brauchen in § 5 das gleiche ‚Integral, aber mit den Grenzen k, 
und oo und schreiben dementsprechend 
L= 
Mit der ‚Abkürzung ¢ aus (22f) haben wir nach (32b) 
2 kt 


Das Korrektionsglied erster Ordnung wird daher beim horizontalen 
Dipol nach (29), (30) und (32c) 


=— 4) — sing 4). 


(83) 


Dabei wurde gesetzt 
(34) = |k, (- <5<0, vgl. Anm. 1’von 8. 9). 


2h 
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Um das Entsprechende fiir den vertikalen Dipol zu machen, miissen 
wir den in (22c) vorkommenden Ausdruck + + 1 etwas umformen, * 


zunächst ohne Vernachlässigung, indem wir auf gleichen Nenner 
bringen, den Nenner rational machen und die Abkürzung benutzen 
k, 2 
(34a) ur 
Man findet so: 

kf? 
— 

M 2k,* k,? 
(34b) WV +1l= 


Wenn man im Zähler des zweiten Bruches nur das zweite und 
dritte Glied 'beibehält, die von der Größenordnung k,? bzw. k,* sind 
und daher das erste überwiegen, und im Nenner des ersten Bruches 
k,* gegen k,* streicht, so erhält man mit (32) und der zu (22d) 
analogen Substitution 24 =k, u für den Integranden von K in (22c) 


Man multipliziert nun aus, dividiert durch u? — u,? und zerlegt den 
Rest in Partialbriiche. Wenn man dabei nur Glieder der Ord- 
nung k,/k, und k,*/k,? beibehält, ergibt sich: 


-tu__2kı* k, 1 
(34c) et Fars +u + du. 


Dabei ist, die Zweideutigkeit von u, so entschieden worden, daß 


kı 
V- kr! 


gesetzt wurde, wobei /—k,* wie früher mit positivem Realteil ge- 
rechnet wird. Auf Grund "dieser Festlegung erhält der eigentlich 
noch hinzukommende Partialbruch mit dem Nenner u+ u, einen 
verschwindend kleinen Faktor. Außer diesem Gliede kann für das 
zunächst Folgende auch das erste Glied in der Klammer von (34c) 
gestrichen werden, weil es von der zweiten Ordnung in k,/k, ist. 
Das zweite Glied dieser Klammer führt auf das frühere Integral L 
in (32a) zurück und unterscheidet sich von ihm nur formal durch die 
Wahl der Integrationsvariabeln (früher A, jetzt u). Das dritte Glied 
derselben Klammer endlich liefert das Integral 


(34d) u=+ 


(84e) Gm fen 
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Zusammenfassend erhalten wir so aus (34c) als Korrektion 
erster Näherung beim vertikalen Dipol: 


(35) Re = Re {7 — 2Re 


Das hier eingeführte Integral G läßt sich auf ER Funktionen 
zurückführen. Wir setzen zu dem Zwecke 


(u - Up) = 
Dann wird zunächst ohne Vernachlässigung 
(35a) G=e-iw e-t 
— 


Wegen der Kleinheit von u,, Gl. (34d), können wir hierfür ohne 
numerischen Fehler schreiben (vgl. aber Zusatz 2) 


(35b) : G= fe 
Eine führt auf 


(36) G= + fier + fe + fe. 
-if 
Das zweite und vierte Glied zusammengenommen sind gleich der 
mit negativem Vorzeichen genommenen Eulerschen Konstanten C. 
Das dritte Glied ist gleich 
— log(— ig) = 
Das erste Glied wird mit —t=is 


ö 0 


Im ganzen haben wir also statt (36): 


3 
d d in 
(86a) 


-- Cig) + 
Hier bedeuten wie üblich: 
= f sin s (Integralsinus) , 
(86 b) 
Ci(t) = C + loge — f (1 — cos s) (Integralkosinus) . 
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Aus (85) und (36a) folgt somit für unser Korrektionsglied in 
definitiver Form!) 


Re K} = Re = (Cie + >) 


= Re {#1} += Tal ma 
Wegen der Bedeutung von d vgl. (34). 


(37) 


§ 5. Die Erdabsorption bei Dipolantennen. 


Den Teil der Strahlung, der durch die Ebene z=h—e im 
Zeitmittel emittiert und von der Erde absorbiert wird, bezeichnen 
wir wie in der Einleitung mit S,,. = — S_ ini Abuse ihn beim 
horizontalen Dipol nach Gl.(24). Die dort angedeutete eckige Klammer 
hat folgenden Inhalt: 


Dabei wurde in der letzten Zeile zu deren Vereinfachung schon von 
der Differentialgleichung der Besselschen Funktionen Gebrauch ge- 
macht. Wir können aber auch die erste Zeile auf die entsprechende 
Form, nämlich ayf das Produkt J(Ar)J{lr) reduzieren, sofern wir 
bier die in (24) vorgesehene Integration nach r ausführen. Wir 


betrachten also das Integral f Xrdr, wo 
0 


(38) X= JI" (ind (nals 


Aus den Differentialgleichungen für J(Ar) und J (Ir) folgt durch ge- 
eignete Multiplikation derselben: 


co po 
(88a) | Xrdr=Al] (lr rdr— and’ dnjar. 


Da hier das zweite Integral rechts verschwindet, hat man: 


f Xrar =A Lf JanJ (rar. 
0 0 


1) Beim vertikalen Dipol sind die Rechnungen bei Niessen soweit 
geführt, daß man die Ergebnisse für das Korrektionsglied erster Ordnung mit 
den unseren vergleichen kann. Es besteht keine volle Übereinstimmung, was, 
wie im Zusatz 3 genauer diskutiert werden soll, auf eine Vernachlässigung bei 
Herrn Niessen zurückzuführen ist. 


23 
! 
| 
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Hiernach gelingt die Zurückführung des dreifachen Integrals in. 
Gl. (24) auf ein einfaches!) genau wie früher in GI. (20). 

. Man erhält, wenn man wieder, was bequem ist, unter dem 
Zeichen Re die Vertauschung von +i mit — i vornimmt: 


Wenn man für f,...f, die Ausdrücke (16), ay einführt, aus- 
multipliziert und beachtet, daß ~ 


(40) m=—m", füra<k, und By -+nfür2>k, 


so kann man zunächst ohne Vernachlässigung das Ergebnis wie folgt 
zusammenfassen: 


ronal 
k, 


(41) 


ky kı 
+ | | ®@ada, 
+f fe, 


1) In den dem gleichen Zweck dienenden Rechnungen von Herrn Niessen | 
fehlt der Beitrag, der von f Xr dr herrührt (Aun. d. Phys. 82, vgl. insbesondere 


8.452 u. 453 daselbst), Diesem Beitrag entépricht in seiner Bezeichnung der 
Ausdruck wo 


I J” An) I’ 


oo 
i, = f 
QO. 


(Dabei haben wir fiw sein r ‘unser A und für sein ¢ unser I beibehalten). 
Nach G. N. Watson, A treatise on the theory of Bessel functions, Cam- 
bridge 1922 Gl. (1) S. 405 ist 


Nun schließt Niessen so 


1 


” ı>1, 


also i,Al+i,= 0 sowohl für A< als auch für A>1. Der Fehler liegt 
offenbar daran, daß die Betrachtung für 1 = nicht mehr richtig ist, weil die 
zu differenzierende Funktion i, an dieser ‘Stelle eine unstetige — hat. 
Annalen der Physik. 5. Folge. 41. 
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@, (a) = Re | 


Der Gang der Rechnung wird in Zusatz 4 gebracht. 
Wir wollen die drei Teile der Reihe nach mit I, II, III be- 


zeichnen. 
I ergibt sich nach (30) und (32c) ohne weiteres 


* 700 2k,* i 
(41b) I=—Refilf} = Bel 
wo Ö die gleiche Bedeutung hat, wie in (34). 
Bei der Ausrechnung von II machen wir die Vernach- 


lässigungen (31), (32) und & 4|<|u,|; wir erhalten 


= Vk,?— |k,| sin 3, 


| | — My |? 1 

| N | ky |* (kt — 2%) 
Es ist dann konsequent auch A? in der eckigen Klammer von ®, zu 
streichen. Damit wird 


(41d) Il = + k,* sind. 


SchlieBlich erhalten wir nach (32) und (34) einfach 


(41 ¢) | 


hy 
de) Ter J add = 


4 
Als gesamte rdabecnytion ergibt sich somit 
Sie verschwindet, wie wir rts # = — > oo (vgl. die Einleitung) 
und wird bei Annäherung der Antenne an die Erde wie £-? un- 


endlich. 

Andererseits ist die Erdabsorption veim vertikalen Dipol durch (20) 
gegeben. Hier kann man ohne weiteres nach Gl. (3) und (3a) zu einem 
einfachen Integral übergehen. Man erhält mit Rücksicht auf (9a) 


ck, A* 
M\ 
| N} wy 


(43) 


1).sind ist negativ, vgl. Anm. 1 von 8.9. 
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Dies wird, wenn man ausmultipliziert und (40) beachtet: 
ky 3 

M aal M |* di 

N “mf 
Der erste Teil, den wir mit IV bezeichnen wollen, hängt zusammen 
mit dem Integral K in (22c), wenn wir in letzterem die untere Grenze 


durch k, ersetzen, was wir wieder durch die Schreibweise an- 
deuten wollen. Wir haben dann 


(438) 


(44) IV = 2Re {iKy} 
= 


wobei zweite Glied rechts daher rührt, daß K mit dem 


Faktor = M +1, IV mit dem Faktor + definiert ist. Dieses 


zweite Gliea verschwindet aber, weil das hinter Re{i stehende 
Integral reell ist. Im ersten Gliede können wir die in Gl. (35) durch- 
geführte Umformung von K verwerten, sofern wir nur rechts L 
durch Ly aus Gl. (32c) ersetzen und in G, Gl. (34e), die untere 
Grenze —i mit 0 vertauschen und dementsprechend statt G 
schreiben G,. 

Wir erhalten so 

k, {sind AN 

Das Integral G,” läßt sich nach dem Verfahren von S. 15 auf die 
Form (36) bringen, wenn man dort nur die untere Grenze —i¢. 
durch — Zu, ersetzt und dementsprechend in (36a) — iu, statt ¢ 
schreibt. Es ist, aber für hinreichend kleines Zu, nach (36b) 


Si(—7ifu,)=0, 
Ci(— = C— + loge + log u, 
— F + log —id 
vgl. auch (34d) Aus (36a) folgt daraufhin 
(44b) =—C— log th 
Indem wir dies in _ einsetzen, erhalten wir 


Wir gehen jetzt auf (43a) zurück. Den zweiten Teil dieses 


Ausdrucks wollen wir mit V bezeichnen und wie folgt auswerten: 


(44c) IV =— 4k,? 


4 
i 
f 
ay 
\ 
ER 
: 
| 
4 . 
\ 


20 Annalen der- Physik. 5. Folge. Band 41. 1942 


Nach den Definitionen (8a) und mit der Vernachlässigung (32) kann 
man schreiben 


M 2k,? 


und hieraus ohne weitere Vernachlässigung mit aa Substitution (22d): 


Mi k, “sind 
(454) 


Damit wird 


- d)udu 


Der Integrand ist, wenn man RB: und im Zähler ( Glieder . 
mit (Er er) vernachlässigt 
u— 20 + 


Hieraus durch Integration zwischen 0 und 1: 


— -20+ — a? log la? + 4] 


Setzt man hier die Bedeutung von a und b aus (45a) ein und be- 
achtet, daß in (45b) vor dem Integral der Faktor k, /|k,| steht, so 
kann man in (450) das zweite und dritte Glied streichen. Das 


vierte gibt — log—*- , das letzte tg 3(F +) 
Somit schließlich ‘ 


(dd V=- 4k,? a (2 + log — ($+ 2) te a) 


Im ganzen haben wir fiir den vertikalen Sender nach (43), (44c) 
und (45d) 


28 k, 
(46) . (Fr — log + 4) sin 
mit der 
(46a) 4 + - (F +4) 


Ebenso wie is horizontalen Sender, Gl. (42), wird auch beim 
vertikalen Sender die Erdabsorption für <—> 0 unendlich wie £2; 
für k, —> oo verschwindet sie natürlich auch jetzt.- 


| (45 c) 


x 
k, 
i 
4 = 
k 
| b = cos d. 
| 
1 |k,| (u — a)* + b* h 
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Wir machen noch darauf aufmerksam, daß (46) für ¢—-> 00 ein 
logarithmisches Unendlichwerden der Erdabsorption anzeigt, welches 
natürlich nicht real sein kann. Wir werden dies in Zusatz 2 
aufklären. 

§ 6. Diskussion und graphische Darstellung der Ergebnisse 


Wir schicken den Grenzfall unendlich gut leitender Erde voraus. 
In den dimensionslosen Variabeln 


25 
28 
haben wir nach (23) und (29) für k, —» oo 


- cost), 


(47) = = 424, 


(48) 
— + cos). 


Für £—» oo ergibt sich 


(48 a) sy = 


Dies entspricht der Ausstrahlung eines Hertzschen Dipols in den 

unbegrenzten Raum, die natürlich unabhängig von der Achsen- 

richtung des Dipols, also die gleiche für sy und sy ist. Dabei 

möge man bemerken, daß der Faktor k,* im Nenner von Gl. (47) 

das Rayleighsche A*-Gesetz widerspiegelt. 
Für haben wir 4 


Bei der Annäherung des Dipols an die Erdoberfläche verdoppelt 
sich also die Ausstrahlung des vertikalen Dipols gegenüber dem 
Werte ?/, bei weggedachter Erde. Dagegen verschwindet die Aus- 
strahlüng des horizontalen Dipols. Man kann dies ganz schematisch 
erläutern durch Konstruktion der mit entgegengesetztem Vorzeichen 
genommenen Spiegelbilder der Einzelpole, vgl. Abb. 2a und 2b: In 
Abb. 2a heben sich die inneren Pole auf und es bleibt ein Dipol 
von doppeltem Hebelarm übrig. In Abb. 2b heben sich die beiden 
Polpaare auf. 
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In Abb. 3 tragen wir den Inhalt der Gl. (48) auf. ¢ = x ent- 
spricht nach Gl, (47) dem Abstande h = as t= 2r dem Abstande 


h= 3 usw. Die Schlängelungen in beiden Kurven deuten die 


Abb. 2. Schema für die Annäherung des Dipols an die vollkommen leitende 
Erde nach der Spiegelungsmethode: a) beim vertikalen Dipol verdoppelt sich 
der Hebelarm, b) beim horizontalen Dipol verschwindet er 


Abb. 3. Die dimensionslosen Strahlungscharakteristiken s, und 8, in Ab- 


hängigkeit von der Lage des Dipols h bzw. von {= . bei vollkommen 
leitender Erde, nach (48) berechnet 


Interferenz der Strahlung von Antenne und Spiegelbild an. Sie 
nehmen mit wachsendem £ ab, wobei sich beide Kurven der Ge- 
raden s=?/,, Gl.(48a), anschmiegen. Der Abstand zweier aufeinander- 
folgender Maxima oder Minima beträgt ungefähr 2x in der Skala ¢, 


; 
de 
vo 
be 
(4 
46 
b) be 
un 
(49 
der 
4 dan 
3 fäh 
| 
der 
Ko 
wir 
das 
(60) 
4 Fer 
Wir 
Laer. 
Der 
ein 


Sommerfeld u. Renner. Strahlungsenergie und Erdabsorption usw. 


das ist A/2 in der Skala h gemessen, wie es nach der Spiegelungs- 


vorstellung sein muß. 
Wir gehen zum Korrektionsglied erster Ordnung über, zunächst 


beim vertikalen Dipol. Nach (23), (37) und (33) haben wir 
4s, = + (sin(¢ — 6) — fcos( — ö)) 
+ sin — cosd (si 3) 


Für ¢-->o wird 
(49a) Asy = 0, 


(49) 


da Ci(oo) = 0.und Si(co) = > wird. Die unendlich entfernte Erde 


beeinflußt also bei beliebiger Leitfähigkeit die Ausstrahlung nicht. 

‘Fir ¢=0 wird nach (36b) Si(0)=0 und Ci()=C + logg 

und daher nach (49) bei konsequenter Näherung: 
(49b) Asy = + — C= log ¢) sind] + F 0s 3]. 
Der endliche .Wert 

sy= + 
den wir in (48b) bei £ =O für k, -> oo ausgerechnet und in Abb. 3 
dargestellt hatten, ist also irreführend. Wenn man die endliche Leit- 
fähigkeit berücksichtigt, kommt eine mit abnehmendem £ unendlich 
anwachsende Korrektion hinzu. 

In Abb. 4 zeichnen wir die Kurve s, aus Abb. 3 noch einmal 
hin, ergänzen sie aber durch den an */, anschließenden Abschnitt 
der positiren Ordinatenachse, die den Limes für ¢—+» 0 des 
Korrektionsgliedes bedeute. Dieses Korrektionsglied berechnen 
wir beispielsweise für Seewasser. Hier ist in der Formel (1) für k,? 
das ¢-Glied bei nicht zu großem w überwiegend. k,* kann also als 
rein imaginär angesehen werden; daher wird in (34) 


. 1 
(50) 


Ferner wird für 4 = 40 m, vgl. den Eingang von § 4, 
k, 


Wir haben also nach (49b) 
Asy = + oe 
Der in Abb. 3 gezeichneten bei */, endigenden Kurve ist hiernach 


ein sehr steil ansteigender Ast zu überlagern, der sich bei ab- 
nehmendem ¢ eng an die Ordinatenachse anschmiegt. Für größere L£, 


gis : 
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z.B. {=.r, ist die Korrektion unbeträchtlich. Das Genauere ist 
aus Abb. 4 abzulesen. 
In dieser ist überdies eine Kurve für 
AK 
4’ 10 
eingetragen. Sie zeigt, wie selbstverständlich, ein weniger enges 
Anschmiegen an die Ordinatenachse für {—» 0 und für gréBere ¢ 
eine stärkere Abweichung von dem Grenzfall des vollkommenen Leiters. 


Abb. 4. Die Strahlungscharakteristik s, in Abhängigkeit von h. Kurve a aus 
Abb. 3 übertragen und durch die Ordinatenachse oberhalb */, ergänzt, ‘Kurve ! 


7 nach (48) und (49) berechnet 


und ¢ für (721 = 100 und 10 mit 8=— 
Es würde keine Schwierigkeit machen, unsere Kurven für 
andere Bodenarten, z. B. für Süßwasser, zu ergänzen, wo wegen des 
erheblich kleineren o das s-Glied in Gl. (1) bei gleichen Wellen- 
längen nicht zu vernachlässigen wäre. Jedoch wollen wir unsere 
Darstellung nicht mit weiteren diesbezüglichen Einzelheiten belasten. 
Beim horizontalen Dipol ist das Korrektionsglied erster Ord- 
nung 48, nach (29) und (33) gegeben durch 


(61) A sy Fink -9- Leong - 9) 
Wir haben also wie vorher für ¢—» co 
(51a) 48,=0 
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und fir &—>0 
1 1 


Das Korrektionsglied wird sonach auch hier unendlich; der frühere 
endliche Wert bei vollkommenem Leiter, nämlich sz = 0 ist inkorrekt. 
In Abb. 5 ist die Kurve für s, aus Abb. 3 durch die positive 

Sy 


Abb. 5. Die Strahlungscharakteristik s,, in Abhängigkeit von A. Kurve a aus 
Abb. 3 übertragen und durch die positive Ordinatenachse ergänzt, Kurve b 


und ¢ für “+l — 100 und 10 mit d=—— nach (48) und (51) berechnet 


Ordinatenachse — und für d= — durch die beiden Kurven 


hei ‘alae überraschende Unendlichwerden von 4s für 
£—>0 rührt, wie wir jetzt zeigen wollen, von demjenigen Teil der 
Gesamtstrahlung her, der der Erde zugestrahlt und von dieser 
absorbiert wird, also von.dem im vorigen Paragraphen berechneten 
Betrage S,,,. Wir definieren s,,, aus S,,,, in derselben Weise wie 
in (47) s aus S und überzeugen uns, daß in der Differenz As — ars. 
sowohl beim vertikaler wie beim horizontalen Sender die für <—>» 0 
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unendlich werdenden Glieder 1/{? und log£ herausfallen. Wir 
haben nämlich nach (49b) und (46) bzw. nach (51b) und (42): 
1 . 

- 4) |sin 3| + 50082], 

(As — Sabs. =(. 

Im Zusatz 5 wird auf einem mehr grundsätzlichen Wege gezeigt 
werden, daß die in den Halbraum oberhalb der Antenne ausgestrahlte 
Energie (S, in unserer früheren Bezeichnung) in jedem Falle, also 
auch für <—>0 endlich ist, so daß die Differenz zwischen der 
Gesamtausstrahlung S und der Erdabsorption S,,. sicher endlich 
bleibt, nicht nur in den mit k,/|k,| proportionalen Gliedern erster 
Ordnung, sondern in beliebiger Ordnung. Wir können diese Differenz 
etwa als „Nutzstrahlung“ bezeichnen. 

Für die der Antenne zuzuführende Leistung (vgl. den folgenden 
Paragraphen) kommt es natürlich nicht auf diese Nutzstrahlung, 
sondern auf die Gesamtstrahlung einschließlich der in der Erde 
absorbierten Energie an, so daß für diese unsere Kurven in Abb. 4 
und 5 maßgebend bleiben. 

Es könnte so scheinen, als ob wir die Erdabsorption statt, wie 
wir es getan haben, durch Integration der Strahlung über die Ebene 
z= hs vielmehr durch Integration über die Erdoberfläche z = 0 
berechnen müßten. Demgegenüber weisen wir darauf hin, daß beide 
Integrationen genau denselben Wert ergeben müssen, was man schon 
daraus erkennt, daß sich zwischen den Ebenen z=h—e und z=0 
kein absorbierendes Material befindet und daß die Ausstrahlung 
durch den unendlich fernen Zylindermantel von der Höhe k für alle 
endlichen £ und daher auch im Limes £ —> oo verschwindet. 
Rechnerisch überzeugt man sich davon folgendermaßen: Wenn man 
in (43) nicht die in Gl.(9a) für z = h — s angegebenen Werte von f, 
und f, einsetzt, sondern die folgenden für z= 0 geltenden Ausdrücke 

(er +e ng), | 
so kommt man zu demselben Ausdruck wie in (42a). Dies gilt für 
den vertikalen Dipol. Entsprechend kann man beim horizontalen 
Dipol verfahren und sich überzeugen, daß auch hier der EnergieduB 
durch die Ebene z= 0 der gleiche ist, wie der durch die Ebene 
zeh-e 

§ 7. Vergleich des idealen Dipols mit einer realen Antenne. 
Der Strahlungswiderstand 

Wenn wir von der Ohmschen Wärme im Antennendraht und 

"etwaigen Funkenverlusten absehen, wird die der Antenne zugeführte 
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Leistung L vollständig in Strahlung umgesetzt, die ihrerseits besteht 
aus der dem Luftraum übermittelten Nutzstrahlung und der in der 
Erde absorbierten Strahlungsenergie. Wir setzen dementsprechend 
in den Gl. (23) und (29) 
S=L 

and erhalten mit der in Abb. 4 und 5 dargestellten Bedeutung 
von Sy und sy 


(53) L=<ck,* A*s,, L = 
2 Vv 


Diese Gleichungen lehren, wie die in (6) und (11) eingeführten 
„Dipolstärken“ A "nd B zu wählen sind, um sie einer Vertikal- 
oder Horizontalantenne gegebener Leistung L anzupassen. Ins- 
besondere seien A, und B, diejenigen Dipolstärken, die der Leistung 
von 1kW = 10 aes enteprechen. Machen wir noch k = —— 


4 in Zentimetern gemessen, so haben wir 

und bei beliebiger in Kilowatt ausgedrückter Leistung: 
(53b) A=A, VL, B=B,YyL. 
Mit den so erhaltenen numerischen Werten von A und B liefern 
die Gl. (9) und (10) die Feldkomponenten in gewöhnlichen Gauss- 
schen Einheiten, also H in Oerstedt und E in elektrostatischen 
Einheiten oder nach Multiplikation mit 300 in Volt/cm. 

Es ist in der Praxis üblich, von der Strahlung S überzugehen 
zu dem Strahlungswiderstand R mittels der Beziehung 
(54) S=RP=+RJ,, 
wo J, die Stromamplitude bedeutet. Hier muß, wenn R in den 
gewöhnlichen elektromagnetischen CGS-Einheiten berechnet werden 
soll, J, in denselben Einheiten ausgedrückt werden. Da wir den 
Ohmschen Widerstand in der Antenne vernachlässigt haben, setzen 
- wir auch 
(54a) L=+R4J,’. 
Während die A, B nach (53a,b) nur von 4 und L und wegen sy, sy 
noch von ¢ und k,/k, abhängen, hängt R durch Vermittlung von J, 
auch noch von der Antennenlänge / und von der Stromverteilung 
ab. Als solche brauchen wir nur die gleichférmige Verteilung über 
die ganze Lange der Linearantenne zu betrachten, also den Fall 
starker Endkapazitäten, da sich jede andere Verteilung in bekannter 
Weise durch einen „Formfaktor“ auf jene zurückführen läßt. 
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Um den Zusammenhang zwischen unseren bisherigen Pol- 
‚stärken A, B und der in (54), (54a) gemeinten Stromamplitude J, 
festzustellen, erinnern wir daran, daß Hertz in seiner klassischen 
Arbeit!) vom Jahre 1888 die unseren A, B entsprechende Größe 
als „elektrisches Moment des Dipols“ bezeichnet. Schreiben wir 
für letzteres ea, so ist e die hin- und herschwingende, bei Hertz 
punktförmig gedachte Ladung, a ihre Schwingungsamplitude. e wird 
bei Hertz elektrostatisch gemessen. Betrachten wir irgendeinen 
zur Dipolachse senkrechten Querschnitt, so können wir, im Interesse 
der Folgenden, e als diejenige Ladung definieren, die diesen Quer- 
schnitt in der Zeit r/2 einer halben Schwingungsdauer einmal 
passiert. 

In unserem Fall ist a gleich der halben Antennenlänge, also 
gleich 1/2. Bezeichnen wir den Antennenstrom mit 


J= J,sinot = J,sin 


so ist die Ladung, die einen Querschnitt in der Zeit r/2 durchsetzt 


t 2/2 
t 
(55) farms, ata 


Dieser Strom bzw. seine Amplitude J, sollte aber elektromagnetisch 
gemessen werden. Um ‘ihn ins elektrostatische Maß zu übertragen, 
müssen wir ihn mit c multiplizieren. Das Hertzsche e ist also 
zu setzen ast i 

Das Hertzsche elektrische Moment ea wird also 


Daher wird auch 
(55a) 
Wenn wir die Ausdrücke für L aus (53) und (54a) gleichsetzen und 
für A, B den Wert (55a) eintragen, so hebt sich 4.44 beiderseits - 


fort und es bleibt wegen k, = as: 


ae | 84 
als Wert des Strahlungswiderstandes in elektromagnetischen CGS- 
Einheiten. Wir gehen schließlich durch Multiplikation mit 10° zu 


1) Die Kräfte elektrischer Schwingungen behandelt nach der Maxwell- 
schen Theorie, Ges. Werke Bd. II, S. 151. 
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Ohm über und fügen den obengenannten Formfaktor & hinzu, der 
bei unserer gleichförmigen Stromverteilung den Wert 1 hat (bei sinus- 


förmiger Verteilung ist a = 2). Dann wird 


(66a) R = 1202" (+) { ro 


8y 
Für eine weit von der Erde entfernte Antenne ergibt sich daraus, 
mit sy = 2, Gl. (48a) 


(56b) R = 802? (5) 2. 


Denselben Wert gibt im Anschluß an die genannte Arbeit von Hertz 
z.B. Zenneck?) an für eine im freien Raum befindliche Antenne. 

Wesentlich ist aber, daß unsere Gl. (56a), über diesen Sonderfall 
hinausgehend, durch Vermittlung der Faktoren s,, sz den Einfluß 
der Antennenrichtung und ihres Abstandes von der Erde, sowie den 
Einfluß der Bodenbeschaffenheit enthält. Besonders das starke An- 
wachsen von sy und Ss, bei kleinen Antennenhöhen h, vgl. Abb. 4 
und 5, kommt, wie wir aus (56a) ersehen, in vollem Betrage auch 
im Strahlungswiderstand R zur Geltung. Diese Abbildungen können 
direkt auch als Darstellung des Strahlungswiderstandes R angesehen 


werden, wenn man nur an der Ordinatenskala den Faktor 120 x* (+) ‘ 


aus Gl. (56a) anbringt, der von Antennenlage und Bodenbeschaffen- 
heit unabhängig ist. 

Man versteht in diesem Zusammenhange auch, daß der Wirkungs- 
grad einer nahe über der Erde befindlichen Antenne günstig be- 
einflußt werden kann, wenn unmittelbar unter der Antenne, wo das 
Feld und daher auch die Erdabsorption am stärksten sein würde, ein 
Netz metallischer Leiter angeordnet wird, das den effektiven Wert 
von k,/|k,! und daher auch von $,,.. heruntersetzt. 

In der ausgebreiteten Lehrbuchliteratur, die im wesentlichen 
mit der Hertz-Zenneckschen Formel (56b) übereinstimmt, ist’ von 
dieser Beeinflussung des Strahlungswiderstandes durch den Erdboden 
und durch die Lage der Antenne gegen denselben nicht die Rede. 
Wir sehen daher in der Klarstellung dieser Umstände ein Haupt- 
ergebnis unserer Arbeit. 


Zusatz 1. Formeln für das Medium II 


Wir haben in § 1 die Potentiale und die daraus abgeleiteten Feldstärken 
nur für das Medium I (Luft) angegeben. Wir wollen hier die entsprechenden 
Formeln für das Medium II zusammenstellen: 


1) J. Zenneck , Lehrbuch der drahtlosen Telegraphie, Stuttgart 1916, 8.49. 
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Zunächst die Verallgemeinerung der GI. (4) für ein Medium der Wellen- 
zahl k, : 
> > 
| 4 =0, 
| § rot If. 
Beim vertikalen Dipol spezialisiert sich, wie wir wissen, der Vektor i auf die 
eine Komponente JZ,. Ihre Darstellung im Medium II ist: 


N 


co 
0 


Der Ausdruck wird einfacher wie der entsprechende im Medium I, und zwar 
deshalb, weil im Medium II keine primäre Erregung vorhanden ist und weil 
man nicht zwischen z < h und z > h unterscheiden muß. 
i 

(6) H,=- 
Man kann nun leicht verifizieren, daß die aus (8) bzw. (8) nach (5) bzw. (5’) 
zu erhaltenden Komponenten E, und H, für z = 0 stetig ineinander übergehen. 

Beim horizontalen Dipol braucht man, wie früher zwei Komponenten JZ, 
und IZ, des Hertzschen Vektors. Diese sind 

Te 
My + Mr’ 
Pdi 

+ 


B 

_ i- 


Nach (4’) wird 
(10°) 


(2 
th, \r 99 
10 1 

, ar 

Man überzeugt sich daraufhin auch hier ohne Schwierigkeit, daß die vier Rand- 
bedingungen für z = 0 erfüllt sind. (In E, und H, enthalten alle Glieder den 


gemeinsamen Faktor cosg, in ZH, und H, den gemeinsamen Faktor sin p). 


Zusatz 2. Verschärfung der Näherungen beim vertikalen Dipol 
für sehr große Höhen 
In Gl. (35a) wurde wegen der Kleinheit von u, die untere Grenze —{(i + u,) 
durch — i¢ und ferner exp. (— £w,) durch 1 ersetzt. Die erste Vernachlässigung 
ist immer zulässig, die zweite RE nur dann, wenn |{w|<1 ist. In 


unserem Zahlenbeispiel = ist dus für > 54 nicht mehr der Fall. 


0° IT, 0° IT, 
= * 4 
E, = cosy + + k,? cos IZ. , 
E sin @ I. sing 1 
N 
= 
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Wir wollen daher hier noch das Ergebnis für Re{iK/k,*} angeben, das man 
erhält, wenn man die Vernachlässigung |{w,| nicht macht. 
Man hat dann statt (37) 


wo @ wie früher durch (36a)-gegeben ist. 

In der gleichen Weise wollen wir auch die Berechnung von Re {ixe} 
ergänzen, die wir in Gl. (44) zur Bestimmung der Erdabsorption brauchten. 
Es wird: 


Dabei bedeuten Lr bzw. G> wie früher die Ausdrücke (32c) bzw. (44b). Das 
noch stehengebliebene Integral erledigt sich durch Reihenentwicklung von.e~ ¥ 
oder bei großem {u, durch Zurückführung auf den Logarithmus. 

Man sieht nun unmittelbar, daß das logarithmische Unendlichwerden 
von Sy, beim vertikalen Dipol, Gl. (46) für ——> ©, auf das wir dort auf- 
merksam machten, nur scheinbar ist. Da nämlich der Realteil von u, positiv 
. ist, verschwindet exp. (— {w,) für {—> ® so stark, daB'der in G,® vor- 
kommende log { unschädlich gemacht wird. Übrigens bemerke man, daß dieser 
Logarithmus überdies durch das zusätzliche Integral in der ( ) der letzten 
Gleichung kompensiert wird. 


Zusatz 3. Vaugpateis mit Niessen, Ann.d. Phys. [5] 22. 8. 162 
und 24, 8. 31. 1935 


Niessen gibt für rein imaginäres k,* und |{w|<1 als Korrektionsglied 
1. Ordnung beim vertikalen Dipol, in unsere Bezeichnungen umgeschrieben, 
an (a. a. O. 8. 180): 


k, 

(ciw + 816 

Das stimmt mit unserem Ergebnis in (37) und (33) überein, wenn wir 
dort 3 =— 7 (entsprechend einem rein imaginären k,*) setzen bis auf das 


Glied — 1 in Fe eckigen Klammer, das bei uns fehlt. Wir können im folgenden 
zeigen, daß dieser Unterschied durch eine Vernachlässigung bei Niessen zu- 


stande kommt, die für h < si sicher unzulässig ist. 
Bei Niessen tritt 8. 174 das folgende Integral auf: 


_wenn wir für sein p unser k,, für sein P unser |k,| setzen. Mit der Sub-. 
stitution Ya? —1 = + yu entsteht: 


U=r 
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u*—2u +2 


{c) U=r udu, x=hryk. 


v 
Hier wird nun bei Niessen öu?/2 gegen 1 vernachlässigt. Daß dies nicht 
immer zulässig ist, soll durch das in unseren Abbildungen 3 und 4 benutzte 
Zahlenbeispiel d = 10~* belegt werden. Nehmen wir u = 100, dann ist das bei 
Niessen vernachlässigte Glied ungefähr gleich 1. Die Umgebung der Stelle 


u = 100 gibt aber für h < a noch einen wesentlichen Beitrag zum Integral, 


weil hier exp. (— xu) keineswegs klein ist. Darüber hinaus können wir direkt 
zeigen, daß diese Vernachlässigung tatsächlich zu dem angegebenen Unter- 
schied zwischen dem Ergebnis von Niessen und unsefem führt. 

Wir gehen auf die Form (b) zurück, machen zunächst die Substitution 
y = Vx* — 1 und spalten U in zwei Teile: 


oo oo 
f ry +9 
In dem zweiten Integral macht man sodann die schon genannte Niessensche 


Substitution y = Zu und erhält 


(e) verfe wo (1 —d + ud). 


2u + 


Wenn man in der aa Klammer — 6 und wö gegen- ı vernachlässigt (5 ist 
nach Voraussetzung etwa 10”* und daher ud etwa 10”? für u = 100), so unter- 
scheidet sich der so aus (e) erhaltene Näherungswert von dem aus (c) durch 
Vernachlässigung von a erhaltenen Niessenschen Näherungswert durch 
das erste Integral in (e). Dieses beträgt 


V2 
Wern man nun das Niessensche Resultat (a) um diesen Betrag korrigiert, so 
wird gerade das fragliche Glied — 1 in der eckigen Klammer von (a) kompensiert 
und die Abweichung von unserem: Ergebnis in (32) und (33) getilgt. 

Die gleichen Bemerkungen gelten auch fiir die in der Uberschrift zitierte 
zweite Arbeit von Niessen, in der die Beschränkung auf rein imaginäres k,* 
fallen gelassen wird. 

Wir sind auf diesen Punkt ausfihrlich eingegangen, weil er fiir den 
Verlauf unserer Abbildungen 4 und 5 entscheidend ist. Nach Niessen würde 
das Glied mit 1/{* in unserer Gl. (49b) wegfallen und fur das logarithmische 
Glied bestehen bleiben. Das würde aber bedeuten, daß das Anwachsen des 
Strahlungswiderstandes mit Annäherung an die Erdoberfläche praktisch be: 
deutungslos sein würde, womit ein Hauptresultat unserer Untersuchung ver- 
loren ginge. 

bb brigens weist Herr Niessen gelegentlich selbst auf Schwierigkeiten 
beim Grenzübergang h —> 0 (divergente Integrale) hin, vgl. a. a. O. Aun. d. 
Phys. 22, S, 176. 
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Zusatz 4. Zwischenrechnung betreffend die Erdabsorption 
beim horizontalen Dipol 


Wir wollen zeigen, wie man von Gl. (39) zu Gl. (41) und (41a) kommt: 
Wenn man in (39) für f, ... /, die Ausdrücke (16), (17) einsetzt und den in (26) 
bereits berechneten Teil von dort übernimmt, erhält man 


ck, Bi 
wo VI—IX folgende Abkürzungen bedeuten 


© 
2 

| 


(a) = VI + VII + VIII + IX + 


4 
ck, Bt 


ky 
1 Pr 
VII | Vi "Pa, 


u + Be 


2 2_ 322 


IX = 2Re ms 


0 


co 
M 


VIII = Re 2h (ur + | Mı — Me 2 1? 
0 


Wenn man zunächst VI und VII'zusammenfaßt, erhält man mit der Be- 
zeichnung ®, aus (41a) ziemlich direkt 


ky 
0 


Für VIII schreiben wir, indem wir den in (41a) angegebenen ersten Ausdruck 
von ®, benutzen: 


V 
Weiterhin empfiehlt es sich, IX in zwei Teile aufzuspalten 
IX=U+V. 


U bzw. V seien diejenigen Integralsummen, die entstehen, wenn wir in der 
Definitionsgleichung von IX die obere bzw. untere Grenze durch k, ersetzen. 
Dieselben Integralsummen treten aber auch in dem Ausdruck von S, Gl, (26) 
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und (26a) auf, die wir bei entsprechender Aufspaltung und Benutzung von (40) 


schreiben können 
48 


—— — 3 

ki 
Infolgedessen fassen sich die beiden letzten Glieder von (a) folgendermaßen 
zusammen ai 


(a) 
Durch Addition von (b), (c) und (d) entsteht nach . 


(f) 2V- fo na-ılama. 


-U+V. 


4,3 
-zh +2V. 


Dies ist aber mit (41) identisch, wenn wir uns noch davon überzeugen, daß V 
nach Gl. (27a) geschrieben werden kann: 


| 


ky 
Zur Vervollständigung der an (41) anschließenden Rechnungen müssen 
wir nur noch zeigen, wie man zu dem an zweiter Stelle in (41a) angegebenen 
Näherungsausdruck von ®, kommt. Dazu genügt es, anzunähern 
Ha Us 


durch —1+ 
Kı + 
und infolgedessen 


-1 durch - 


= Re {— . 
+ Hg ) { 


|? 


Zusatz 5. Über die Nutzstrahlung im Grenzfalle {—> 0 


Wir wollen hier die in $6 aufgestellte Behauptung beweisen, daß S, 
stets und in Strenge endlich ist und es auch im Grenzfalle t——> 0 bleibt. 
Wir gehen aus von dem Ausdruck für den vertikalen Dipol, Gl. (20), (20a): 


48 ] 


Wir setzen für f, und f, die Werte aus (9a) (mit dem oberen Vorzeichen) ein, 
multiplizieren aus und berücksichtigen die ee Es entsteht: 


ki 
V 


M Adi 

N 

Die beiden ersten Integrale sind von h unabhingig. Beim letzten Integral 

kann man zum Limes A—> 0 übergehen, weil die Integrationsgrenzen 0 


und k, sind. Also behält auch dieses für ¢—-» 0 einen endlichen Wert. (Bei 
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: §_ kam das gleiche Integral vor, die Integrationsgrenzen waren aber dort k, 
und ©, so daß man dort unter dem Integral nicht zu {—> 0 übergehen 
konnte) Die analoge, aber etwas umfangreichere Betrachtung beim 'hori- 
zontalen Dipol wollen wir hier unterdrücken. 


Zusatz 6. a) Die zweite Näherung beim horizontalen Dipol 


Bei der Ausrechnung der gesamten ausgestrahlten Energie in $ 3 hatten 
wir das Korrektionsglied mit Re {i L/k,*} bezeichnet, wo L durch (30) und das 
dort vorkommende Fi) durch (27b) definiert waren. Indem man in FA) 
Zähler und Nenner mit M multipliziert und Glieder von der as k,3/k,* 


und k,*/k,* wegläßt, erhält man 
ke 


In §4 haben wir nur das erste Glied berücksichtigt. Wir wollen jetzt den 
Beitrag der zweiten Näherung ausrechnen. Das in Betracht kommende Integral 
kam in § 3 Gl. (28) schon vor, dort'aber mit den Grenzen 0 bis k,, während 
hier durch den Faktor %k,?/k,? die Realitätsverhältnisse geändert sind und die 
Grenzen durch 0 und © ersetzt werden müssen. Man rechnet leicht nach, daß 


SE Die zweite Näherung des ggg: wird daraufhin: 

k,? 
— [@?— 1) sin — 28) + — 20)]. 
LA 


Sante wollen wir die zweite Näherung zu den drei Gliedern I, II, III, die in 
dem Ausdruck (41) für die Erdabsorption beim horizontalen Dipol vorkamen, 
angeben. 

Zu I. Wenn wir für F(4) den Ausdruck dieses Zusatzes nehmen, erhalten 
wir für des Glied zweiter Näherung: 


—y Re +il- i} = 


Zu II. Man überzeugt sich leicht, daß man die Gl. (41e) ersetzen muß 
durch die we. um eine Potenz von k,/|k,| genaueren Werte: 


— 2 2k — 

Damit erhält man als pe 2. Ordnung in II 
ThE J 5 


Zu III. Bei ®, DER sich in Vervollständigung der am Ende von Zusatz 4 
angedeuteten Rechnung als zusätzliches Glied zweiter Näherung 


k’sin?d. 


3* 
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Das liefert integriert als entsprechende Korrektion zu III: 


16 k,? 

15 
Durch Vergleich mit den Gliedern 1. Ordnung sieht man, daß die bei II 
und III auftretenden Glieder 2. Ordnung durchweg zu vernachlässigen sind. 
Das bei I neu hinzukommende Glied spielt erst für ganz kleine Antennen- 
höhen eine Rolle, wenn nämlich 


k,’sin?öd. 


z1 


wird. In unserem Zahlenbeispiel kann man sich daher für h > 5 ohne 
weiteres mit der Rechnung erster Näherung begnügen. 


Zusatz 6. b) Zweite Näherung beim vertikalen Dipol 


In $4 wurde für das Korrektionsglied Re {i K/k,®} in dem Ausdruck für 
die gesamte ausgestrahlte Energie S nur das Glied von der Ordnung A,/]|k,| 
ausgerechnet. Wir wollen jetzt das Glied mit k,?/|k,|* angeben. Es kommt 


von dem in (34e) vernachlässigten Term k,//— k,* her und wird 


Auch bei der Berechnung der Erdabsorption in $5 müssen wir die ent- 
sprechende Ergänzung machen. In Gl. (44) kommt der Ausdruck 2 Re {iKy} 
vor, den wir dort ebenfalls nur in erster Näherung ausgerechnet haben. Das 
zugehörige Glied 2. Ordnung heißt: 

Aus dem in $5 mit V bezeichneten Gliede entsteht kein Glied 2. Ordnung. 

Durch Vergleich mit den entsprechenden Gliedern erster Näherung sieht 
man, daß hier beim vertikalen Dipol die Glieder 2. Ordnung überhaupt keine 
Rolle spielen. 

Wir möchten nicht verfehlen, auch an dieser Stelle den Herren 
W.Moser und W.Magnus bestens zu danken für wertvolle An- 
regungen und Diskussionen. 


München, im November 1941. 


ki ps sin 20. 


(Eingegangen 29. November 1941) 
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Réntgenographische Fouriersynthese 
von metallischem Magnesium 


Von R. Brill, C. Hermann und Cl, Peters 
(Forschungslaboratorium Oppau der I. G. Farbenindustrie A.-G.) 
(Mit 3 Abbildungen) 


In Fortsetzung unserer Arbeiten über den Zusammenhang 
zwischen Elektronenverteilung und chemischer Bindung?) geben wir 
nachstehend eine ausführliche Beschreibung unserer Versuche und 
der erhaltenen Ergebnisse bei metallischem Magnesium’). 

Es ist selbstverständlich, daß bei einem Metall eine zwei- 
dimensionale Fouriersynthese nicht zum Ziele führt, weil alle Metall- 
gitter irgendwelche dichteren Packungen darstellen, so daß man keine 
Projektionsrichtung finden kann, bei der nicht Überlappungen der 
Elektronenwolken störend wirken und das eigentliche Ziel unserer 
Arbeit, nämlich eine klare Erkenntnis über die Verteilung der Elek- 
tronen zwischen den Metallatomen oder -ionen zu erhalten, un- 
möglich machen. Es mußte daher eine dreidimensionale Fourier- 
synthese vorgenommen werden, was zur Folge hat, daß die Durch- 
messung einer Zone nicht genügt, sondern sämtliche überhaupt 
auftretenden meßbaren Reflexe heranzuziehen sind. 

Als Untersuchungsmaterial dienten uns gut ausgebildete Magne- 
siumeinkristalle, die aus der Dampfphase gewonnen waren (Abb. 1). 
Bevorzugt trat die Basis auf, die jedoch bei den meisten Kristallen 
Krümmungen aufwies. Diese waren uns bereits ein Hinweis auf 
Unregelmäßigkeiten in der Sekundärstruktur der Kristalle, so daß zu 
erwarten war, daß sie eine für unsere Zwecke hinreichende Mosaik- 
struktur aufweisen. 

Ein geeignetes, aus dem Präparat Abb. 1 ausgesuchtes Kristall- 
individuum wurde beiderseitig parallel zur Basis angeschliffen, so daß 
eine Kristallplatte von gleichmäßiger Dicke (0,90 mm) entstand?). 
Diese Platte wurde in geeigneter Weise auf dem a.a. O. beschriebenen 
Goniometertisch befestigt. Der Tisch erlaubte zwei Drehungen um 


1) R.Brill, .G.Grimm, C. Hermann u. Cl. Peters [5] 34 S. 393. 
1939, in Hinkunft mit a. a. O. zitiert. 

2). Vorläufige Mitteilung vgl. Naturw. 26. 8. 479. 1938. 

'3) Für das Schleifen dieser Kristallplatte sind wir Herrn Prof. Dr. Himmel, 
Heidelberg, sehr zu Dank verpflichtet. 
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eine horizontale und eine vertikale Achse (vgl. Abb. 1, a. a. O.), so 
daß ein justierter Kristall in beliebige, definierte Richtungen eingestellt 
werden kann. Zur Messung der Intensitäten wurde die a. a. O. be- 
schriebene Apparatur verwendet. 

Da die intensitivsten Basisreflexe (insbesondere 0002) einen 
Intensitätswert ergaben, der mit dem theoretischen, aus der Hartree- 
kurve entnommenen sehr gut übereinstimmte, durfte unsere Kristall- 


Abb. 1. Magnesiumkristalle 


platte als extinktionsfrei betrachtet werden. Ein weiterer Hinweis 
auf die Freiheit von primärer oder sekundärer Extinktion war uns 
die ziemlich große Breite der Reflexe. (Nach unseren Erfahrungen 
kommt es nur sehr selten vor, daß ein Kristall von vornherein 
extinktionsfrei ist.) 

Die Form des Kristalles brachte es mit sich, daß fast alle 
Interferenzen entweder als Durchgangsreflexe oder als Oberflächen- 
reflexe gemessen werden konnten. Nur solche Reflexe, bei denen der 
eintretende oder austretende Strahl die Basis nahe tangential streift, 
ließen sich mit unserer experimentellen Anordnung nicht erfassen 
(Schattenwirkung der Haltevorrichtung) und ermöglichten auch keine 
definierte Berücksichtigung der Absorptionswirkung. Sie wurden 
daher nicht gemessen. Ihr Wert ließ sich aus dem experimentell 
bestimmten Verlauf der F-Kurve interpolieren. Bis zu einem 
Reflexionswinkel von 30° wurden die Reflexe mit an Steinsalz mono- 
chromatisierter Strahlung gemessen. Die Intensitäten von Inter- 
ferenzen mit größerem Beugungswinkel wurden ohne Monochroma- 
tisierung mit der natürlichen Strahlung der Röntgenröhre vermessen, 
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wobei jedoch solche Reflexe, bei denen eine Koinzidenz mit 6-Strahlung 
möglich war, von vornherein ausgeschlossen wurden. 

Der Wert des Absorptionskoeffizienten wurde experimentell an 
der für unsere Messungen verwendeten Kristallplatte .zu u = 6,59 
‘bestimmt), Die gemessenen Intensitätswerte und die mit diesem u 
berechneten F-Werte sind in der Tab. 1 zusammengestellt. Dabei 
wurden folgende Formeln benutzt: 


1. Für Durchgangsreflexe (Lauefall) 
J=Q-1-cscO, 
2. Für Oberflächenreflexe (Braggfall): 
J =Qlesc@, Ste. 


Dabei ist J = ae die integrierte Intensität des Reflexes. 


Ferner bedeuten: 
l = Dicke der: Kristallplatte. 

0 =» + der Eintrittswinkel. 

0, = «— # für Durchgangs-, + — « für Oberflächenreflexe, der 
Austrittswinkel (d.h. der Winkel zwischen dem betreffenden 
Strahl und der Kristalloberfläche). 

= Glanzwinkel. 
a = Neigung der Netzebene gegen die Kristalloberfläche. 
F = Strukturfaktor je Atom Mg. 


= Zul(csc, + cse 9,). 


1 
é= (csc O, csc @,). 
1 1 + cos*23 - F? das Reflexionsvermögen der 


sin 24 


untersuchten Netzebene pro cm des Primärstrahlweges ?). 


1) An einem polykristallinen Magnesiumblech wurde bestimmt u = 7,01, 
während der Wert nach Jönsson 7,62 beträgt. 
2) Für monochromatisierte Strahlung müßte diese Formel etwas geändert 


werden: An Stelle des Polarisationsfaktors + (1 + cos* 29) sollte hier stehen 
1 + 23 cos*2 9, 
1 + cos*2%, 


wo & = 7°15’ der Glanzwinkel für Mo Ka an (200) von Steinsalz ist. Wie 
man sich leicht überzeugt, bedeutet diese Vernachlässigung bei 9 = 20° einen 
Fehler von weniger als 1°/,, bei 30° von 2°/,, was noch in den Bereich der 
Streuung der Intensitätsmessungen fällt. 
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v = Atomvolumen von Mg. 

e,m = Ladung und Masse des Elektrons. 
c = Lichtgeschwindigkeit. 
1 = Röntgenwellenlänge. 

Die Tab. 1 läßt erkennen, daß alle Reflexe, die von hohen 
Ordnungen der Basis herrühren, auf höhere F-Werte führten als 
entsprechende Prismenreflexe. Da die Basisreflexe, wie erwähnt, 
stets als Oberflächenreflexe, die Prismenreflexe jedoch als Durch- 
gangsreflexe gemessen werden mußten, schien uns der Unterschied 
zwischen den Intensitäten dieser beiden Reflexarten mit der Art der 
Messung zusammenzuhängen. In der Tat ergab sich, daß auch an 
der Basis verbotene Reflexe, wie z. B. 00011 oder 00013, mit nicht 
monochromatisierter Strahlung gemessen werden konnten. Da solche 
Reflexe mit monochromatischer Strahlung nicht auftreten, müssen 
sie durch das Spektrum des kontinuierlichen Untergrundes verursacht 
werden. Dann sollte aber ein analoger Effekt auch für die Prismen- 
reflexe erwartet werden. Da er nicht auftritt, nehmen wir an, daß 
er bei den Basisreflexen durch den langwelligen Anteil des Kontinuums 
verursacht wird, der bei den Durchgangsreflexen unverhältnismäßig 
viel stärker geschwächt wird als bei Oberflächenreflexion. Wir haben 
daher die höheren Ordnungen der Basisreflexe für unsere Auswertung 
nicht verwendet, sondern statt dessen eine glatte Kurve durch die 
für die Prismenreflexe bei größerem Winkel gefundenen Intensitäts- 
F-Werte gezogen. 

Die Fouriersynthese wurde in zweifacher Weise durchgeführt. 

1. Mit den auf einer glatten Kurve interpolierten F-Werten 
der Tab. 1b. 

2. Allein mit den mit monochromatischer Strahlung gemessenen!) 
Intensitätswerten unter Anbringung einer „Temperaturkorrektion“ zur 
Unterdrückung des Abbrucheffektes 


Die unter 1. angeführte Fouriersynthese lieferte die in Tab. 2a 

dargestellten Ergebnisse. Der Untergrund ist nicht völlig befriedigend 
und zeigt noch eine geringe Nullpunktsunterschreitung am Punkte E. 
Der Grund dafür ist in dem hier nicht auskorrigierten Abbruchs- 
effekt zu sehen. 
1) D.h. für die Reflexe mit 5 < 1,3. Nur bei dem völlig aus der Reihe 
fallenden Reflex (10 14) wurde der interpolierte statt des gemessenen F-Wertes 
eingesetzt. Hier dürfte der oben erwähnte Fall eines undefinierten Absorptions- 
weges infolge zu kleinen Austrittswinkels aus der Kristallplatte vargelegen 
haben. 


2) Vgl. a. a. O. sowie Ztschr. f. Elektrochem. 46. S. 425. 1940. 5 
y's 
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Tabelle 1a 


Gemessene Intensitäten J, Strukturamplituden F für das Mg-Gitter (bezogen auf 


Foooo = 12); Atomfaktoren Fy, für das Mg-Atom, vergl. mit den auf Zimmer- 


temperatur") umgerechneten Werten F, nach Hartree und F, nach Pauling 


1/d 


‚385 
‚769 
1,154 
1,539 


0 
1 
2 
3 
4 
5 
6 
7 
20 
22 
24 
26 
20 
21 
22 
23 
24 
25 
26 
0 
1 
2 
3 
4 
5 
0 
2 
4 


so 


888 


| 


| f 


| 1808 | . 412, | 


1) Charakteristische Temperatur 9 = 295°. 
2) Intensitäten, die im Braggfall, d.h. als Oberflächenreflexe gemessen 
‘werden, sind durch B gekennzeichnet. 


en > 
ed 0006 1053 B. 
ied 00 6 214 B 8,365 — F, ; 
ler 0001 2,92 B 4,740 | 4,740 8,21 a 3 | 
2,308 ‚112 B 10 4 
ht 14 2,693 B 316 1,03: 2,46 
‚409 41.0 4.230 | ‚06 
680 7,42 106 7.985 8,44 
850 7,68 3,137 6.974 7,99 8,84 
ab 1,027 139 1.38 8,28 
1,209 ‚300 B ‚138 3.476 5,558 7,10 
1,394 0 B 3,080 3.557 4,20 5,84 
1 ‚624 2.090 2,10 ° „27 
die 5,260 | 5,260 606 | 628 
ts- 721 3 = 
‚746 ‚67 2,77 1,70 3,50 oes 
‚817 8,40 y 70 5.540 1,70 fier 
92 2,14 422 | 5, 5,23 
‚924 23 ‚110 55 
1,054 3,16 ‚370 4,740 5,02 ‚ol 
| © 360 | 30 
n') 958 0 2,205 | Ex = 
1 
1,35 ‚160 ‚427 ‚62 cl 
E 0 1,082 | | 268 
0 1108 170 3,058 
9 1,249 - : 0 2,68 Bes: 
ihe 04 | "885 | 1.976 
ns- 5 0 3 ‚295 | ‚750 1 
gen | 
1 
59 
—— 
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Tabelle 1a (Fortsetzung) 
hikt | 17a | 7.18 F > 
3360 1,874 ‚045 ‚548 ‚48 
6060 | 2,161 075 219, |  ,219 25 
Tabelle 1b. Interpolierte Werte Fy, und ¥ (Grundlagen für die Rechnung 1) 
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42 
; 
hikt | 
0002 385 
4 ‚269 
6 1,154 
Se 8 | 1,539 
10 | 1,924 
ae 12 | 2,308 
14 | 2,693 
3 as 1010 ‚362 | 
4 "850 ‚12 
5 | 1,027 03 
6 | 1,209 
s | | 147 3080 | 1,082 | 302 | 302 
9 | i768 | 2 | 1,148 | 2,67 | 2,67 
10 | 1957 | 751 
11 | 2,146 „0 61 
12 | 2336 | ‚14 
Ber 13 2,527 *05 05 
14 | 2,717 01 
BR - 1120 | ,624 | 6,14 1,98 
a 2 ‚233 | 5,25 1,49 
4 991 | 3,50 90 
ae. 6 | 1,312 | 1,96 "46 
a 8 1,661 ‚98 ‚17 
Re 10 | 2,022 „42 02 
12 | 2,301.) 11 
14 | 2,765 0 1,00 
1,69 
2030 | 5,34 90 
u; 1 ‚46 | 5,15 1,39 
ee 2 3817 | 4,63 ‚68 
3 924 | 3,91 ‚93 
4 1,054 | 3,16 2 
5 1,202 | 2,39 2,07 
6 1,361 1,80 90 ‚22 
Be: 7 | 1,528 | 1,81 | 1,13 ‚25 
= 8 | 1,200 | 90 | ,45 07 
ae 9 1,875 ‚62 ‚54 ‚06 
as 10 | 2,054 | ,39 „20 01 
Be 11 | 2,235 | ,22 ‚19 
= ae 12 | 2,418 09 ‚05 
13 | 2602 | ‚08 | ‚08 
14 | 2,788 0 0 
Be: 


ow 
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hikı | ı/a Fug| F hikl | 1/4 F |1jd F 
4040 |1,442/1,55| ‚8 4260 |1,908| ‚57 | 29 4480 |2,498/ ,06 |,06 
1 | 1,455 | 1,52 | 1,32 1 | 1,917] ,55 | ,48 2 | 2,528] ‚05 | ,05 
2 11,498 |1,41| ‚70 2 | 1,946| ‚52 | ‚26 4 |2,613| ‚03 |,03 
3 | 1,554 | 1,24 | 1,07 3 11 ‚45 | ‚39 6 2,751] ‚01 | ,01 
4 |1,635|1,04| ‚52 4 |2,057| ‚38 | ‚19 
5 1947 ‚84 5 31 a {5350 
6 |1,847| ,66| , 6 | 2,230| 23 | ,1 10 
7 |1,973| ,48| ‚42 7 | 2,335] ,15 | ‚13 
8 12,109| ‚33] ‚17 8 | 2,451) ‚08 | ‚04 2 2'553 ‚04 02 
9 |2,253/ ,21| ‚18 9 |2,575| ‚04 | ‚03 3 | 2'590| 503 | 403 
10 | 2, ‚10| ,05 10| 2,709 ‚01 | ‚01 4 | 2,639 02 01 
1 5 /2,791 | ,01 |,01 
5160 |2,008| ,43'| ,22 6 2,775) 0 | 0 
1 [2,016 | ‚42 | ‚36 
3250 [1,572|1,19| ‚60 2 12,044 ‚40 | 20 6280 |2,600| ‚08 | ,02 
1 |1,584 | 1,16 | 1,00 3 | 2,089] ‚35 | ‚30 1 |2,607 | ‚03 | ,03 
2 |1,618/1,07| ‚54 4 |2,150| ‚30 | ‚15 2 | 2,629| ‚03 |,02 
3 | 1,674|. ,95| ‚82 5 | 2,226) ‚23 | ‚20 3 | 2,663 |. ,02 | ,02 
4 |1,750| ,83| ‚42 6 |2,315| ‚16 | ‚08 4 12,711 | ,01 |,01 
5 |1,842| ‚68| ‚50 7 |2,418| ‚09 | ‚08 5 |2,772| 0 | 0 
6 |1,950| ,51| ‚26 8 | 2,529| ‚05 | ‚03 
7 |2,070| ,37| ,32 9 |2,651| 02 | 02 130 01 | 01 
8 |2199| 18 —— 2 |2)749 | ;o1 |,01 
9 12,339| ‚14| ‚12 6060 |2,161| ,25 | ‚25 
10|2,484 | ‚07| , 2 |2,197| ,21 | ‚21 
11 | 2,636} ,02| ‚02 4 |2,296| ‚17 | ,17 
6 | 2,452] ‚08 | ‚08 
4150 1,652 1,00 | 1,00 8 | 2,655} ,02 | ‚02 
2/1 91| ‚9 
4 |1,823| ;70| 4370 |2,193| ,26 | ‚13 
6 |2,015| 2| , 1 |23,201 | ‚25 | ‚22 
8 |2,258| ‚20 2 |2,226 | ‚23 | ,12 
10|2,535| ,05| ‚05 3 12,268| ‚20 | ‚17 
4 |2,325| ‚15 | ‚08 
= 5 | 2,395| ,10 | ,09 
5050 |1,803| ,74| ‚37 6 |2,479| ‚07 | ‚04 
1 |1,813| ‚72| ‚62 7 12,574) ‚04 | ‚03 
2 |1,843| ,67| ‚34 8 |2,679| ‚02 | ,01 
4 |1, 25 5270 |2,267| ,21 | ,21 
5 [2,043] ‚40| ,35 2 |2,284| ‚18 | ‚18 
6 12,140) ,30) ‚1 4 |2,380| ,11 | ,11 
y D D 8 727| ‚01 1 
9 2,500 = 2, 
6170 !2,365| ‚12 | ‚06 
1 | 2,372] ‚12 | ‚10 
3360 |1,874| ,62| ,62 2 12,395| ,10 |, 
2 !1,913| ‚57| ,57 3 |2,434| ‚09 |, 
4 |2,025| ,41| ,41 4 |2,486| ‚06 | ‚03 
6 |2,200| ‚25| , 5 12,553 | ‚04 | ‚08 
8 |2,424| ‚09| ‚09 6 | 2,631! ‚02 | ‚01 
10 |2,685 | | ‘01 7 |2721| ‚oı | 
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Die Fouriersynthese nach 2. führt zu Ergebnissen, die in Tab. 2b 
niedergelegt und auszugsweise bereits an anderen Orten ausführlich 
dargestellt und diskutiert wurden’). 


Abb. 2. Gemessene F-Werte von Magnesium, verglichen mit der glatten Inter- 
polationskurve (ausgezogen) sowie mit den theoretischen Kurven nach Hartree für 
Mg-Atome und Mg* * -Ionen (gestrichelt) und nach Pauling für Mg-Atome (strich- 
punktiert). Die Basisreflexe, die nicht monochromatisch gemessen wurden und den 
im Text erwähnten Intensitätseffekt zeigten, sind als ausgefüllte Kreise dargestellt 


Hier soll noch eingegangen werden auf die Berechnung des 
Kugelradius, der gerade die 10 Elektronen des Mg* *-Ions einschließt. 
Die Anzahl N der Elektronen innerhalb einer Kugel vom Radius r 
in einer kugelsymmetrischen Ladungsverteilung N (r) ist gegeben 
durch das Integral: 


N(r) = 42] o(a)2*dz. 
J 


Für unsere Rechnung, deren Ergebnisse noch Meß-, Abbruchs- 
und Abrundungsfehler enthalten können, dürfte die Näherung ge- 


1) Vorläufige Mitteilung vgl. Naturw. 26. S. 479. 1938. 
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Tabelle 2a. Rechnung 1 
Dichten in Elektronen. pro AS 


Kantenlänge des dargestellten Dreiecks —“- = 1,85 A 
Höhe 0-c=0Ä 
‚05 25 113 4,93 
5. ,23 


v3 


—,08 
—,17 E; H; I; K; L; 0% 
Höhe !/,c = 0,650 4 


*) Die Buchstaben bezeichnen die entsprechenden Punkte in Abb. 3. 


| 
2b 
ch 
05 6,20 A; B: 1 
N C; D;G 
06 ,10 ,17 24 1,04 410 26,34 = 
‚3 ,08 ‚14 ,20 39 174 855 33,20 
131 . 4,08 
10 13 52° 179 4,88 
8 07 7: 27 88 2,10 E 
10 13 ,33 1,06 
MM: 
‚il ,13 ‚20 ‚24 
‚08 ‚08 ‚12 ‚21 
01 07 09 13 — 
01 10 ‚07 
04 ,10 07 | 
‚07 2 
o1 02 — 
ir 9. ,09 : „13,22 - 1,03 
8 ,i0 ‘20° - 6- - 1,12 
‚08 1 09 ‚09 0 ‚15 2 23 23 8 82 - 
It 10° 08- 15°: 28° - 28: - 45 
: ‚06 ‚07 2 ‚14 ‚18 
1 13° ‚14 
‚09 10° 09 ‚09 ‚09 3 
10: 1 
9 - 10 ‚10 
1 
,09 ’ 
‚0,0 
‚12 
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Tabelle 2a (Fortsetzung) 
Höhe */,,¢ = 0,975 A 


Höhe */,,¢ = 0,217 A 


3,25 18,03 36,14 
u 
18,00 


46 
= ,08 ‚08 ‚08 ‚09 13 ‚19 y21 
‚08 ‚08 ‚08 ‚09 ‚09 ‚15 ‚18 ‚24 
‚08 ‚07 ‚09 ‚08 ‚10 ‚16 ‚22 
En: ‚06 ‚06 ‚09 ‚07 ‚12 ‚19 
,06 ‚04 ‚09 ‚09 ‚07 ‚16 „21 
‚04 ‚05 ‚09 ‚08 ‚10 ‚19 
5 8 10 ,06 ,15 
‚08 ‚08 ‚09 ‚0 
03 00 
04 09 
07 
oo 16 1,00 
eo ,09 ‚12 ‚16 ‚22 ‚19 
‚09 yll ‚15 yal 335 
‚10 ‚16 ‚22 
07 08 ‚10 ‚21 
3" 10’ 16 
‚03 ‚09 ‚12 
06 „10 
02 
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Tabelle 2a (Fortsetzung) 
Höhe !/,,c = 0,433 A 


‚10 ‚08 ‚il ‚24 ‚36 1,58 3,77 6,35 
‚09 ‚09 ‚14 ‚18 386 2,36 5,23 
‚08 ‚18 ‚36 1,46 3,30 = 
‚09 ‚10 ‚12 ‚14 ‚22 ‚2 1,99 3,85 oe © 
‚08 ‚12 ‚10 „el 30 1,14 2,37 
‚10 Ad ‚12 ‚23 ‚48 1,46 
‚08 ‚08 ‚il ‚17 ‚26 ‚68 1,58 Be 
‚07 ‚09 jll ‚22 ‚28 ‚85 
„07 ‚10 ‚14 ‚22 „34 
,06 ‚10 ‚10 17 
‚07 ‚10 ‚19 ‚19 4 
‚09 „10 ‚10 „17 
‚09 ‚09 ‚07 ‚12 ‚18 
07 ‚06 ‚08 ‚14 = ae 
‚07 ‚07 ‚10 
‚06 ‚03 ‚08 ‚12 
‚07 ‚08 yll 
‚07 ‚08 4 
07 
07,08 : 
‚11 
‚13 | 
: Höhe */,c = 1,300 A — 
‚12 ‚10 ‚05 ‚09 ‚09 ‚08 ‚10 ‚12 
‚09 „10 ‚09 ‚09 „0 
‚07 ‚10 ‚10 ‚08 ‚10 ‚10 ‚07 
‚09 „il 07 ‚07 ‚il ‚09 
‚10 ‚12 ‚08 ‚08 ‚12 ‚10 
ll ‚10 ‚07 ‚10 ‚il 
‚08 ‚08 ‚08 ‚09 ‚11 
‚10 ‚08 ‚07 ‚08 ‚10 
‚07 ‚07 ‚07 ‚07 
,06 ‚05 ‚06 ‚05 ‚06 aa 
‚03 ‚05 ‚05 ‚03 ue 
‚03 ‚05 ‚03 
,01 ‚04 ‚04 ‚ol 
‚02 ‚04 ‚02 
‚03 ‚03 
‚03 ‚03 ‚03 
‚03 ‚03 > 
‚04 A 
‚05 ‚05 2 
‚06 
2 
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Tabelle 2b 
Rechnung 2 
Höhe O-c=04 


08 E;H, I; K; L; O 


Höhe */,,¢ = 0,325 Ä 


1,53. 3,67 
91 


48 

,18 ,15 ,12 ,71 214 528 9,17 -11,01 A; B; 
N 13 17 Al 122 3,44. 7,18 . 10,36 C; D; 
15 ,13 25 202 497 863 
16 ,13 ,17 39 110 306 636 9,17 
14,14 25 60 1,71 4,14 
4 18 35 243 5 

_ ’ ’ ’ 
i Be: 14 15 24 50 130 309 5,32 - 
15 21 39 91 2,06 
14 16 25 48 1,12 .2,19 
‚13 29  ,60 1,5 
13 19 34 - 
‚10 gie ‚23 ‚42 
15 4 

09 10 "25 
1 15 
08 08 10 
08 08 
‘ 

‚15 ‚12 1,45 3,45 - 5,99 

15” 15" 20° 216. 3,68 
AB r 0 1,23 2,42 
‚15 ‚18 31 ,68 1,45 . 
a ‚14 ‚16 ‚24 ‚46 91 
‚52 
13,14 ‚20 

‚12 
10 ‚10 08 ‚09 08 
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Tabelle 2b (Fortsetzung) 
Höhe */, c = 0,650 A 


DP 


3388 

888 8 8 

33 


® 
SB 
8838 


8 
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ll ‚09 ‚08 13 ‚30 65 1,34 22 4 
‚10 ‚09 ‚09 18 42° 93° 2,65 ‘2 
10 09” ’ 1,26 2,08 
11 ’ 18 ,40 ‚84 1,56 220 
10%» 12’ 26° 4° 1,07’ 1,74 a 
‚1 ‚10 16 34 69 : 
’ ‚10 ‚12 ‚22 41 ‚65 
,13 ‚10 15 27 46 
jll ‚13 ‚20 31 
Al 213 ‚el 
12’ 09" 
07 
A 
’ 
08 „23 ‚sg ‚63 ‚34 = 
,09 20 ‚28 47 ‚69 
‚10 ‚18 ,30 ‚43 y 
,09 
1 ‚15 19 ‚25 ‚4 
28 
10 16 
‚08 yll ; : 9 
11 
’ 
m 
‚06 
‚07 
07 ‚07 
4 Y 
| xum 


8 
1838388 
ES 
28335555 
833835355888 
SE 83 85558 
835383 88835 
38888 

33288 

8 


Höhe */,,¢ = 0,975 A 


Tabelle 2b (Fortsetzung) 
Höhe */,¢ = 1,300 A 


3 


50 
‚5 
5,48 
1° 42 
= ‚4 - 
7 ‚38 
7 
0 8 9 10 ,09 
eee ce ie 08 0 ‚10 0 ‚10 09 ‚08 09 10 10 0 10 
109 m 09 ‚10 99 Pi ‚0 
06° ;06 06 
05" 05" 05 
a 0,4 ‚04 
05 
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nügen, bei der dieses Integral durch eine Summe von g-Werten in 
äquidistanten Punkten n-d ersetzt wird: 


sofern die Punktdistanz d nicht zu groß gewählt wird. 


Abb. 3. Magnesiumgitter mit Einzeichnung der Linien, deren Elektronen- 

dichten in Tab. 3 zusammengestellt sind. Das schraffierte Dreieck oben in 

der Abbildung stellt das Gebiet dar, dessen Dichte in Tab. 2a und 2b als 

„Höhe 0“ angegeben ist. Die übrigen Dreiecke liegen in entsprechender Höhe 
senkrecht über oder auch unter dieser Fläche 


Nimmt man als Grundlage dieser Rechnung die Dichte auf der 
Linie A—P—B mit d =-0,133 A, so ergibt sich: 
Für Rechnung 1: N (6,5 d) = 9,81; N (7,5d) = 10,23 
und durch lineare Interpolation: N (6,91 d) = 10,00. 
Für Rechnung 2: N (8,5 d) = 9,89; N (9,5 d) = 10,49 
und N (8,68 d) = 10,00. 
Die Kugel mit 10 Elektronen hat also 
für Rechnung 1: den Radius r = 6,91d = 0,919 Ä und das Vo- 
lumen 3,25 A’, d.i. 14°/, des Kristallvolumens V, pro Atom, 
das 23,0 A® beträgt; 
für Rechnung 2: den Radius r = 8,68d = 1,154 A und das Vo- 


lumen 6,44 A’, d.i. 28°/, von V,. 
4* 


N 
J 
X 
\h 
j 
| 
| 1 | 
| 4 
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Tabelle 3 

Elektronendichten nach den Rechnungen 1 und 2 längs einigen ausgezeichneten 
Linien im Mg-Gitter. Die Dichten g, (©) sind zugleich die Ordinaten zu den 
in den früheren Mitteilungen veröffentlichten Kurven 


A—R—B C—F—-D H-—-M-N-0 
zÄ & | @ zÄ 


A,B 0 |66,20/11,01 C,D 0 
‚133 |52,71|10,36 1 
‚267 


E88 |? 


= 


23 


“ 


iit 


BBES 


S18 


S88 | 


- 


Da dieses Volumen in der Rechnung 2 durch eine bedeutende 
Rechentemperatur vergréBert ist, hat der Wert, der aus der Rech- 
nung 1 folgt, die gréBere physikalische Bedeutung. Er ist daher 
auch in unseren früheren Mitteilungen publiziert worden. 

Auf eine Wiederholung der übrigen schon früher mitgeteilten 
Ergebnisse verzichten wir und stellen nur zusammenfassend fest, 
daß ein verhältnismäßig glatter Elektronenuntergrund existiert, der 
einer gleichmäßigen Verschmierung von etwa 1,75 Elektronen pro 
Magnesiumatom über das ganze Gitter entspricht. Außerdem ergab 
sich, daß im Gitter ein ausgezeichneter Punkt auftritt, an dem eine 
‚Elektronenanhäufung und in dessen Nachbarschaft Punkte mit 


TG 
736,14) 306) +,07] 
_ 
‚400 | 8,55) 6,36 433| 6,35| — 1186| 04 
‚34| 2,75 4,14 ,850| 1,431 265 1,582] ‚06 
1,31 243 — | 99 1,978 118 
‚s00| 52) 1,30 ‚975 22) 358 2,373] ,45 
‚27 69 1,300} ,07 2,768) 8,56 
1,067 | ,‚23| ‚39 1,625| ,00 B,C 3,164/66,2011 
1201| 316] 1,733; — | ‚08 3,560] 8,56 
1334| 317 1,9501 511 3,955) 50}. 1,34 
1,468; ,11) ,13 2,167) ,10) — 4,350} ‚14| ,23 
1,601| 12] — | 18 4,746] 12] 7 
F 2600, ,18 
4—E—F-G L—E-K 
> 
A 0 |66,20) 11,01 2,705 | ıı LEEO |-,17 
‚225 |33,20| 9,17 2,930 ,14 ‚217 ‚02 
451 | 483| 532| 3,156 | ;06| ‚14 325 | — 
676 | 1,12| 2,19| 3381 | [16 433 | 413 
902 | ‚z2| F 3607 | 305| [18 650 | 712 
4127 | 713| 3832 | 715 367 | — 
1,352 | 307} 4057 | ol 975 | ‚oz 
1578 | 302} 4,283 | 1,300 | {06 | 
E 1903 |-,17| 4508 | 771 
2029 | 4,734 | 1,13] 214 
es 2,254 ‚04| ,08 4,959 | 4,93| 5,28 
2,479 | 5,184 | 33.27| 9,17 
2,705 | ‚ı2| 311|*@ 5,410 | 66,20| 11,01 


35 


2% 
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Elektronenmangel vorhanden sind. Wieweit dieser Effekt reell ist, 
müssen wir dahingestellt sein lassen. Theoretische Gründe!) würden 
nicht gegen das Auftreten eines solchen Punktes sprechen. Außer- 
dem hat bereits Brindley*) aus Messungen der atomaren Streukurve 
von Ruthenium auf Elektronenanhäufungen im Gitter geschlossen. 
Ferner. sind unlängst Untersuchungsergebnisse von Griffith und 
Harris?) an Li mitgeteilt worden, aus denen sich ebenfalls auf der 
Würfelkante des Li-Gitters eine Elektronenanhäufung ergibt. Diese 
unabhängigen Befunde anderer Forscher sprechen für die Realität 


des von uns gefundenen Effektes. 


1) Vgl. A. Sommerfeld, Naturw. 28. S. 769. 1940. 
2) G. W. Brindley, Nature 140. S. 461. 1937. 
3) R. L. Griffith u. P. M. Harris, Phys. Rev. [2] 59. S. 110. 1941. 


Ludwigshafen a. Rhein, 10. November 1941. 


(Eingegangen 15. November 1941) 
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Bemerkungen zum Vorzeichen der transversalen 
thermomagnetischen Effekte im Falle des freien 
Elektronengases. Der transversale Thomsoneffekt 
in anisotropen Metallen und die analogen Er- 
scheinungen im Magnetfeld 


Von Max Kohler 


1. In einer früheren Arbeit!) wurden die transversalen galvano- 
magnetischen und thermomagnetischen Effekte unter der Annahme 
des isotropen und freien Elektronengases betrachtet. Es ergab sich 
dabei zwischen den thermoelektrischen Erscheinungen ohne Magnet- 
feld und dem Ettingshauseneffekt ein Zusammenhang, der experi- 
mentell zumindest qualitativ bei den unter den gemachten Voraus- 
setzungen in Frage kommenden Metallen bestätigt ist. Bei der 
Betrachtung des Ettingshausen-Nernst-Effektes wurde der adia- 
batische Effekt als die dem Experiment zugängliche Größe angesehen. 
Für die adiabatische Ettingshausen-Nernst-Konstante ergibt sich: 

1 

(1) Agra = — = 

Darin ist n die Elektronenkonzentration, ¢, die Fermische Grenz- 
energie am absoluten Nullpunkt und c die Lichtgeschwindigkeit. 
Darnach hätte also Agy, unter den gemachten Voraussetzungen 
stets das negative Vorzeichen. Nun ist aber in einer neuen Arbeit?) ' 
gezeigt worden, daß man bei der Berechnung des adiabatischen 
Ettingshausen-Nernst-Koeffizienten die thermoelektrische Feldstärke 
in den seitlichen Anschlußdrähten berücksichtigen muß. Die Be- 
rücksichtigung dieses Umstandes ergibt Gleichheit des isothermen 
und des so berechneten adiabatischen E.N.-Koeffizienten, sofern nur 
die seitlichen Anschlußdrähte dieselben thermoelektrischen Eigen- 
schaften haben wie die Metallprobe, an der die Messung durch- 
geführt wird. Damit ist aber nicht mehr die Formel (1) maßgebend, 
sondern die Formel: 


(la) Ages = -A 


1) M. Kohler, Ztschr. f. Phys. 118. 8. 37. 1941. In Zukunft als I bé- 
zeichnet. 
2) M. Kohler, Ztschr. f. Phys. (erscheint demnächst). 
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Dies ist die Formel (4i) in I. Darin ist 4’ = (en. wo r 
die von der Energie E abhängige Stoßzeit der Elektronen bedeutet. 
Das Vorzeichen von Agy hängt also vom Vorzeichen von A’ ab. 
Bei den einwertigen Edelmetallen folgt aus dem Vorzeichen der 
Thermokraft ein negatives A’. Daher ist auch Agy bei diesen 
Metallen, in Übereinstimmung mit der Erfahrung, negativ. Bei den 
meisten Alkalien ist A’ auf Grund des gemessenen Thomsonkoeffi- 
zienten positiv anzunehmen. Daher erwarten wir bei diesen Metallen 
ein positives Agy. Messungen liegen nur an Li vor, das auch ein 
positives Agy aufweist‘). Trotz dieser Übereinstimmung fällt aber 
Li außerhalb des Gültigkeitsbereiches unseres Modells. Denn Li 
besitzt eine positive absolute Thermokraft?2), während aus dem 
positiven Vorzeichen von Agy eine negative Thermokraft folgen 
würde. Diese Abweichung ist wohl darauf zurückzuführen, daß bei 
Li Bindungseffekte eine maßgebende Rolle spielen, wie dies auch 
- durch die hohe Ionisierungsspannung des freien Atoms wahrschein- 
lich gemacht wird. 

2. Der gewöhnliche Thomsoneffekt ist in isotropen Körpern so 
definiert, daß in einer vom el. Strom durchflossenen Metallprobe 
bei gleichzeitiger Anwesenheit eines Temperaturgradienten pro 
Volumeinheit und pro Zeiteinheit die Wärmetönung: 

(2) U =—p- (I grad T) 

entsteht. Die Versuchsanordnungen zur Bestimmung des Thomson- 
koeffizienten u sind so getroffen, daß I und grad T parallel gerichtet 
sind. Ist dies nicht der Fall, dann wird nach (2) nur die Kompo- 
nente von grad T in der Stromrichtung für die Wärmeentwicklung 
bestimmend. Anders in anisotropen Medien, z. B. in nichtregulären 
Metaljkristallen außerhalb des Magnetfeldes. In diesem Fall ist die 
Thomsonwärme gegeben durch: 


3 
> oT 
(2a) U =— 
i, 


k=l 
Hat der el. Strom etwa die Richtung der z,-Achse, dann ist die 


gesamte Thomsonwärme nicht mehr durch das Glied — y,, J, ae 


oT 


Ga, und 


gegeben, sondern es kommen noch die Glieder — y,,J, 


1) A. W. Smith, Phys. Rev. 8. 8. 79. 1916. 
- 2) G. Borelius, „Grundlagen des metallischen Zustandes. Physikalische 
Eigenschaften der Metalle“, Handb. d. Metallphysik Bd. I, 1. Teil, 8. 400 
(Leipzig 1935). 
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_ 7 hinzu. Die Konstanten (4 + k) sollen die Koeffi- 


zienten = transversalen Thomsoneffektes genannt werden, in 
Analogie zum transversalen Peltiereffekt und zur transversalen 
Thermokraft. Das Auftreten des transversalen Thomsoneffektes in 
nichtregulären Kristallen hat zur Folge, daß in einem Kristallstab, 
der in der Längsrichtung von einem el. Strom durchflossen ist, 
bei Vorhandensein eines transversalen Temperaturgradienten eine 
Thomsonwärme entsteht. Die Konstanten p,,(i = 1,2,3) sind die 
Koeffizienten des longitudinalen Thomsoneffektes. In isotropen 
Stoffen und in kubisch kristallisierenden Metallen ist u,, = OF + k), 
und alle Bu ( = 1,2,3) sind einander gleich. 

Es sei nun 9 der Winkel des (dem trigonalen, EIERN 
oder hexagonalen Kristallsystem angehörenden) Kristallstabes gegen 
die kristallographische Hauptachse. Der el. Strom sei im rechts- 
händigen rechtwinkligen Koordinatensystem parallel zur z,-Achse. 
Die z,-Achse stehe senkrecht, zur Ebene aus Stab und kristallo- 
graphischer Hauptachse. Dann ist für die vorkommenden Metall- 
kristalle bekanntlich: 

Hy, = COS? + SIN? = = SIN? + COS? 
Pie = = Mes = = 9, “thy = = (Ba — #,) sin cos 

Die Gleichung u,, = 0 besagt, daß ein transversales Temperatur- 
gefälle senkrecht zur Ebene aus Stab und kristallographischer 
Hauptachse keinen transversalen Thomsoneffekt zur Folge hat. 
Die Beziehung #,,= u,, sagt aus, daß man beim transversalen 
Thomsoneffekt Stromrichtung’ und die Richtung des transversalen 
Temperaturgefälles miteinander vertauschen kann, ohne die Größe 
des Effektes zu ändern. Der transversale Thomsoneffekt verschwindet 
für Stäbe parallel und senkrecht zur kristallographischen Maupt- 
achse, denn es ist = (uy — g cos g. 

Zwischen den Komponenten ¢,, und Z7,, der Tensoren der ab- 
soluten Thermokraft und des absoluten Peltiereffektes bestehen 
folgende thermodynamische Relationen: 


m 
Ba) 
(3b) = 


Für die Konstanten u,,( + k) des transversalen Thomsoneffektes 
folgt daraus auch die Beziehung: 


(8c) 
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Bei bekannter Temperaturabhängigkeit von kann daraus 
berechnet werden. In Temperaturen oberhalb der char. Temperatur 
ist &,, näherungsweise proportional zur abs. Temperatur T, daher 
ist dann genähert ¢,, = u... 

Nun zu den entsprechenden Erscheinungen im Magnetfeld. Für 
die Entwicklung der Wärmemenge gilt wieder die Formel (2a). 
Jetzt ist aber zu beachten, daß im allgemeinen %,#y,. Im 
Magnetfeld ist es also allgemein nicht mehr gestattet, Stromrichtung 
und die Richtung des transversalen Temperaturgefälles zu ver- 
tauschen, ohne die Wärmetönung zu ändern. Die Elektronentheorie 
der Metalle liefert bekanntlich auch im Magnetfeld Relationen 
zwischen den Größen u,,($), &,() und J7,,(9). Die Klammer mit 
dem Buchstaben § soll die Abhängigkeit dieser Größen vom Magnet- 
feld § zum Ausdruck bringen. Die Größen ¢,,() und /7,,($) treten 
in folgendem Zusammenhang auf: 

3 


Darin sind F, die FERIEN der el. Feldstärke, (Jy); die Kom- 

ponenten des Wärmestromes, w,, die Komponenten des el. Wider- 

standes und A,, die Komponenten des Wärmeleitfähigkeitstensors. 
Die oben genannten Relationen lauten nun’): 


am? 
) ki A) ) 


Hierin bringt der obere Index (g) oder (u) zum Ausdruck, daß es 
sich um den im Magnetfeld geraden vom ungeraden Anteil dieser 
Größe handelt. Es folgt Besen 


(40) 


(4d) po = — "ar 


Die Beziehung (4c) ist die zu (8c) analoge. Die Größen u” (é + k) 
entsprechen im Magnetfeld dem eigentlichen transversalen Thomson- 
effekt, ähnlich wie die Größen ¢,, (i + k) die eigentliche transversale 


1) M. Kohler, Ann, d. Phys. [5] 40. S. 601. 1941. 
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Thermokraft im Magnetfeld bestimmen. Die Berechtigung zu dieser 
Aussage entnehmen wir einmal der vollkommenen Gleichheit der 
Formeln (3c) und (4c) und aus der in einer früheren Arbeit!) ge- 
fundenen Tatsache, daß die Konstanten = und nM? keinen Beitrag 
zur Entropiebilanz liefern und die ihnen entsprechenden Effekte also 
als reversibel angesehen werden können. 

Die Größen «{), 77% sind nur im Magnetfeld von Null ver- 
schieden, und bestimmen den eigentlichen Ettingshausen-Nernst- und 
Ettingshauseneffekt. Diese Konstanten liefern, wie früher gezeigt, 
einen Beitrag zur Entropiezunahme. Die zugehörigen Effekte sind 
demnach als irreversibel anzusehen. Dasselbe gilt für den durch 
die Konstanten un gegebenen Effekt. Daher auch der erhebliche 


k 
Unterschied im Aussehen der Formeln (4c) und (4d). 
Als Spezialfall betrachten wir den Fall des isotropen Leiters 
im transversalen Magnetfeld. Der el. Strom habe wieder die 
a,-Richtung, das Magnetfeld die Richtung der z,-Achse unseres 
rechtshändigen Koordinatensystems. Dann ist ¢,, ein Maß für den 
Ettingshausen-Nernst-Effekt. Außerdem ist im isotropen Fall: 


ei) =— eS) und ef = ef = 0. Die Ettingshausen - Nernst - Feld- 
stärke F, ist unter Beachtung der üblichen Vorzeichenvorschrift 


gegeben durch: 


Daher ist er =— B,„ H. Nun sei in der z,-Richtung ein Tem- 


peraturgradient angelegt. Dann ist die Wärmetönung 0) 
dx, 


1 


Nach (44) folgt: 


(5) 


oT 
ap (1? + Bex)+ J, 
Man erhält im Magnetfeld bei transversalem Temperaturgefälle also 
eine Wärmetönung, die ihr Vorzeichen ändert bei Umkehr des 
el. Stromes, bei Umkehr des transversalen Temperaturgefilles und 
bei Magnetfeldumkehr. In 1. Näherung ist U proportional zum 
Magnetfeld H. 

In höheren Temperaturen ist nach der Elektronentheorie der 
Metalle Bey näherungsweise temperaturunabhängig, daher folgt dann: 


H-J,- 32. 


1) M. Kohler, Ann. d. Phys. [9 40. S. 601. 1941. 
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Bei bekanntem Ettingshausen-Nernst-Koeffizienten Bgy ist die 
Wärmetönung U in Abhängigkeit vom Magnetfeld zu berechnen. 

In tiefsten Temperaturen, wo nur Gitterunregelmäßigkeiten 
streuen, ist Bey proportional zu T, und damit auch U. 


Zusammenfassung 
Im 1. Teil der Arbeit wird das Vorzeichen des Ettingshausen- 
Nernst-Koeffizienten im Falle des freien Elektronengases näher 
untersucht. Dabei wird die in einer früheren Arbeit angegebene 
Gleichheit des isothermen und adiabatischen Ettingshausen-Nernst- 
Koeffizienten benutzt., Das Vorzeichen des E.N.-Effektes ist abhängig 


von dem Vorzeichen der Größe A’ = (ee. (x = Streuzeit). 


Bei den einwertigen Edelmetallen ergibt sich aus den Messungen 
des Thomsonkoeffizienten ein negativer Wert von 4’, und damit ein 
negativer E.N.-Effekt, in Übereinstimmung mit der Erfahrung. Bei 
den Alkalien ist ein positiver E.N.-Effekt zu erwarten. Das Li 
fällt außerhalb des Gültigkeitsbereiches der angenommenen Näherung. 

Im 2. Teil wird der sogenannte transversale Thomsoneffekt in 
Metallkristallen zunächst bei verschwindendem Magnetfeld betrachtet. 
Darunter wird eine Wärmetönung verstanden, die proportional zur 
el. Stromdichte und zum Temperaturgradienten senkrecht zum Strom 
ist. Ein solcher Effekt kann bei verschwindendem Magnetfeld nur 
in nichtregulären Metallen auftreten. Er verschwindet für Kristall- 
stäbe parallel oder senkrecht zur kristallographischen Hauptachse. 
Der Effekt liefert keinen Beitrag zur Entropiebilanz und ist daher 
reversibel wie die übrigen thermoelektrischen Effekte. Zwischen 
den Koeffizienten der transversalen Thermokraft, denen des trans- 
versalen Peltiereffektes und denen des transversalen Thomsoneffektes 
bestehen die thermodynamischen Relationen (3c). 

Bei Vorhandensein eines Magnetfeldes kann man zunächst die 
in einer früheren Arbeit angegebene Trennung der thermomagneti- 
schen Erscheinungen in die eigentlichen thermoelektrischen Er- 
scheinungen und die eigentlichen verallgemeinerten Ettingshausen- 
oder Ettingshausen-Nernst-Effekte durchführen. Die ersteren sind 
die im Magnetfeld geraden Anteile gewisser Tensorkomponenten, die 
letzteren die im Magnetfeld ungeraden Anteile derselben Größen. 
Die ersteren sind reversibel, die letzteren irreversibel. Für die 
eigentlichen transversalen thermoelektrischen Effekte gelten im 
Magnetfeld genau dieselben thermodynamischen Beziehungen wie 
im Falle verschwindenden Magnetfeldes. Anders dagegen für 
die irreversiblen Effekte. Der dem transversalen Thomsoneffekt 


$ 
! 
=: 
q 
al 
j 
4 
. 
| 
4 


60 Annalen der Physik. 5. Folge. Band 41. 1942 


analoge irreversible Effekt im Magnetfeld besteht in einer Wärme- 
tönung, die proportional zum el. Strom, proportional zum trans- 
versalen Temperaturgradienten und in 1.Näherung proportional zum 
Magnetfeld wird. Diese Wärmetönung ändert also bei Feldumkehr 
ihr Vorzeichen. Der für diesen Effekt maßgebende Koeffizient 
hängt gemäß (4d) mit den Koeffizienten des verallgemeinerten 
E.N.-Effektes zusammen. Der Charakter dieser Relation (4d) ist 
wesentlich verschieden von der entsprechenden Beziehung (4c) zwi- 
schen den analogen reversiblen Effekten. Der Koeffizient dieses 
Effektes ist in höheren Temperaturen temperaturunabhängig, in 
tiefsten Temperaturen, wo nur Gitterunregelmäßigkeiten streuen, 
proportional zu T. 


Berlin, 1. Institut für theoretische Physik der Universität. 


(Eingegangen 8. Dezember 1941) 
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H. Marz. Zur Theorie der Zylinder- und Kugelwellen usw. 


Zur Theorie der Zylinder- und Kugelwellen 
in reibungsfreien Gasen und Flüssigkeiten 


Von Helmut Marx 


(Mit 13 Abbildungen) 


Zusammenfassung 


Untersuchung des mathematischen Verhaltens und der physikalischen 
Bedeutung derjenigen von K. Bechert') angegebenen Lösungen der hydro- 
dynamischen Gleichungen für Zylinder- und Kugelwellen in reibungsfreien 
Gasen und Flüssigkeiten, bei denen die Dichte nur von der Strömungs- 
geschwindigkeit abhängt. — $1. Verlauf der Integralkurven der Differential- 
gleichung (1), I. für Gase, II. für Flüssigkeiten. $2. Beschreibung der den 
Lösungen von (1) entsprechenden Wellenausbreitungen, I. für Gase, II. für 


Flüssigkeiten. 


Die vorliegende Arbeit ist ein Auszug aus der Dissertation des 
Herrn Dr. H. Marx (zur Zeit im Felde); die zum Teil langen mathematischen 
Beweise für die hier vorgebrachten Behauptungen und die Ergebnisse rein 
mathematischer Art der Dissertation erscheinen in den Veröffentlichungen 
des Math. Seminars der Universität Gießen. Die Arbeit schließt an das 
Ergebnis der 8$1—4 meiner Arbeit „Über die Ausbreitung von Zylinder- 
und Kugelwellen in reibungsfreien Gasen und Flüssigkeiten“ ') an, und zwar 
sollen -die dort in §§ 1—4 angegebenen Lösungstypen auf jhr physikalisches 
Verbalten untersucht werden; bei ihnen ist die Dichte g nur von der 
Strömungsgeschwindigkeit « abhängig. Außerdem ist dort vorausgesetzt, daß 

die Zustandsgleichung die Form habe: p =p,+-" ©, (p ist der Druck), 
wobei p,, a, n Konstante sind. Diese (in der Literatur gewöhnlich gemachte) 
Voraussetzung über die Zustandsgleichung bedeutet neben der Annabme g = ¢ (u) 
eine wesentliche, weitere Einschränkung der Lösungsauswahl, denn die Unter- 
suchung der ebenen Wellen in Gasen hat gezeigt”), daß die „Adiabaten- 
konstante“ a? im allgemeinen im Gas von Schicht zu Schicht verschiedene 
Werte hat, also von Raum- und Zeitkoordinaten abhängig ist. 

Zur Feststellung des physikalischen Verhaltens ist die Untersuchung der 
Integral der in der Arbeit!) angegebenen GI. (16) oder (27c) nötig. Die in der 
genannten Arbeit gegebenen allgemeineren Lösungen (vgl. dort § 5 und folgende) 
werden in der vorliegenden Arbeit nicht untersucht. 


Gießen, Institut f. theoret. Physik, 26. Nov. 1941. K. Bechert 


1) K. Bechert, Ann. d. Phys. (5) 89. 8. 169. 1941. 
2) Vgl. K.Bechert; Ann. d. Phys. [5] 89. 8. 370. 1941; ebenda [5] 40. 
8. 207. 1941. 
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In der Arbeit von Prof. Bechert ist auf Grund der hydro- 
dynamischen Grundgleichungen und unter der Voraussetzung, daß 
die Vorgänge reibungsfrej und der Zustandsgleichung p = p, + er 
gehorchend verlaufen, eine gewöhnliche Differentialgleichung 1. Ord- 
nung hergeleitet [vgl. dort Gl. (27c)], deren Lösung die Beschreibung 
der Ausbreitung der Zylinder- und Kugelwellen gestattet, bei denen 
die Dichte eine Funktion der Strömungsgeschwindigkeit u ist. Diese 
Differentialgleichung. lautet: 


dv (k—1) (2w +kin- 1)v 
(1) ‚w* dw kw (w — 1) - 


Hierin ist k = 2 bei Zylinderwellen, k = 3 bei Kugelwellen; n ist 
bei Gasen das Verhältnis c,/c, der spezifischen Wärmen, n=—1 
bei Flüssigkeiten; p, = 0 bei Gasen, p, +0 bei Flüssigkeiten. Die 
Veränderlichen v und w rechnen sich gemäß 


2(k-1) 
—2 3 


2(k—1) 
k 


auf Größen q und g um, aus denen die Beziehungen zwischen 
Strömungsgeschwindigkeit und Dichte 


„[Yke-D 
(1b) u ’ | da | Ay 

Zeit t, und Abstand r von der Symmetrieachse (Zylinderwellen, 
bzw. dem Symmetriemittelpunkt (Kugelwellen) des Vorgangs auf 
Grund von 


(f=): E=yn; 


r n k(n—1) 
ln) 


gewonnen werden können; C und r sind willkürliche Konstanten, 
a ist eine Materialkonstante. 

In § 8 der genannten Arbeit ist die Lösung von (1) für den 
Fall n = — 1 hergeleitet. 


$1. Verlauf der Integralkurven der Differentialgleichung (1) 
Die Differentialgleichung (1) geht mit 


(2) p= 2 ’ z=w—l 
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über in 


(3) 
z w 


Es ist a = 1 für Zylinderwellen, «= 3 für Kugelwellen. 


I. Der Fall der Gase n>1; 6>0) 


Die (w, v}-Ebene sei durch die Geraden v = 0; w= — 8; w=0 
in die Gebiete I, II ..., VIII eingeteilt entsprechend folgender 


Übersicht: 
o 
IlZ Z\iz Z zZ 
w 


Abb. 1. Einteilung der (v, w)-Ebene für Gase. 


a) Zylinderwellen, b) Kugelwellen, 
daher a = 1 daher a = 4 


Dann liest man für den Fall « =+4 (Kugelwellen) aus (3) ab: 


L Wenn v <0, w<P, ist v negativ. 

Il Wenn v<0, —B<w<0O, ist v positiv. 
UI. Wenn v<0, O<w< 1; v negativ. 
IV. Wenn v>0, w>1: v positiv, wenn 


DE =0O wenn 
’ za’ 


. 
negativ, wenn 0 < — 
z 


V. Wenn >0, w<— ß: wv’ positiv, wenn v < =, 


, w* 
v =00, wenn 


v' negativ, wenn v > 
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v>0, —B<w<0: v positiv, wenn >“, 


. w? 
V=0, wenn v=—., 
v’ negativ, 


0<w<1: v positiv, wenn v< =, 
wenn 


v negativ, wenn ® >=. 


v<0, w>1:v positiv. 


v=0: außer won w=0; w=1; 
dort ist v’. unbestimmt. 
w=—f: v=0, außer bei v = 
wo v’ unbestimmt. 
XL w=0: vw =0, außer bei » = 0, wo v’ unbestimmt. 
XI. w=1: v=o; v' unbestimmt bei = 0. 
Daraus ergibt sich qualitativ folgendes Bild für die Schar der 
Integralkurven der Differentialgleichung (1): 


Abb. 2. Die Integralkurven om (1) qualitativ für n > 1 (Gase) 
die gestrichelten Linien stellen v = = vor, sie sind nicht Integralkurven, 


ebensowenig wie w=— 6, w=0. 
Alle anderen ausgezogenen Linien sind Integralkurven 
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Für den Fall « = 1 (Zylinderwellen) ergibt sich in gleicher‘ 
‘Weise ein Kurvenbild, das aus dem obigen durch Spiegelung der 
Ebene w<1 an der Geraden v =0 hervorgeht ( weil nämlich 
für «=1 die Kurve o= = bei w <1 unterhalb der Geraden » = 0 
verläuft). 

Außerdem lassen sich noch folgende Aussagen machen, deren 
Beweis hier zuviel Raum beanspruchen würde: 

Im Gebiet (I,II, III) stellen alle Kurven eine eindeutige 
Beziehung v = v (w) dar und überschreiten dies Gebiet nicht. 
Für w—> —oo hat v den Grenzwert 0; von w=—oo her- 


kommend fallen die Kurven, erreichen bei w=— # ein Minimum, 
steigen an, erreichen bei w = 0 einen Höchstwert, fallen und unter- 
schreiten für w—»1 alle Grenzen; dies gilt für « = = (Kugel- 
wellen); bei & = 1 (Zylinderwellen) ist der Verlauf entsprechend 
spiegelbildlich zur Geraden v = 0. 

Die im Gebiet IV verlaufenden Lösungskurven überschreiten 
dies Gebiet nicht. Jede Integralkurve stellt eine für alle v >0 
erklärte eindeutige Beziehung w = w/(v) dar. Für v— 0 geht w—> 00; 
von ©=0, w=co) -herkommend fällt die Kurve w(v) (w-Achse 
als Ordinate, v-Achse als Abszisse gesehen), trifft auf die Kurve 


v=, hat dort ein Minimum, steigt an und wächst für v->oo 
z 


über alle Grenzen. 
Im Gebiet (7, VI,VII) gibt es vier im Punkte B mit den Koor- 


dinaten w =—f, v= ; : > mindende Integralkurven: Die eine, 
K,, verläuft von B aus bei abnehmendem w bezüglich der v-Werte 
fallend (Kugelwellen) bzw. steigend (Zylinderwellen) und strebt 
für w—» —oo dem Werte v = 0 zu. — Eine weitere, K,, verläuft 
von B aus bei abnehmendem w steigend oberhalb (Kugelwellen) 
® 
bzw. fallend unterhalb (Zylinderwellen) der Kurve = und 
erreicht für w—> —co den Grenzwert |v| oo. — Die dritte, K,, 
verläuft von B aus bei wachsendem w steigend oberhalb (Zyl.) 
bzw. fallend unterhalb (Kug.) der Kurve v= = und unterhalb 
(Zyl.) bzw. oberhalb (Kug.) der Geraden » = 0, und strebt für w—>0 
dem Wert v=0 zu. — Die vierte, K,, verläuft von B aus bei 
wachsendem w steigend, erreicht bei w = 0 ein Maximum, fällt, bis 
sie bei v « < das Steigungsmaß oo hat, bei weiter fallendem v 
z 
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nimmt sodann w ab und für w—>0 strebt v—>0. Dies gilt 
für «= 4 (Kug.); bei «= 1 (Zyl.) ist der Kurvenverlauf ent- 


sprechend spiegelbildlich zur Geraden v = 0. Vgl. Abb. 3, die dem 
Fall der Zylinderwellen entspricht (k = 2, daher « = 1), 

Diese 4 Lösungskurven teilen das Gebiet (V, VI, VII) in 4 Teil- 
gebiete ein, und jede in einem solchen Teilgebiet verlaufende Lösungs- 
kurve verläßt das Teilgebiet nicht [abgesehen von der Möglichkeit 
der Fortsetzung über (w,v) = (0,0) hinaus). Und zwar ist: 

G, das. Teilgebiet, das begrenzt wird durch K,, K, und die 
Gerade v= 0. Alle Integralkurven innerhalb G, verlaufen, bei 
(w, v) = (0,0) beginnend, mit abnehmendem w steigend (Kug.) bzw. 
fallend (Zyl.), erreichen bei w= —£ einen Extremwert und streben 
für w—>—oo dem Wert v= 0 zu. 

Im folgenden wird für die Gebiete G,,G,,G, das Verhalten der 


Kurven beschrieben, wie es fiir a = = (Kugelwellen) gilt. Vgl. auch 


Abb. 2! 

G, wird begrenzt durch K,, K,. Jede Integraikurve stellt eine 
für alle vo >0O erklärte eindeutige Beziehung w = w(v) dar. Von 
(w,v) = (— 00,0) kommend steigt sie bei zunehmendem »v an (v-Achse 
als Abszisse, w-Achse als Ordinate. gesehen), erreicht ein Maximum 


dort, wo v= = ist, fällt und strebt für v—»oo gegen w= —oo. 
z 

G, wird begrenzt durch K,,K,. Von (w,v) = (—o0, 00) kommend 
fällt (bezüglich der v-Werte) jede Kurve in G, bei wachsendem w, 
erreicht bei w=— ein Minimum (min > - steigt an, 
erreicht bei w=0 ein Maximum, fällt, bis sie da, wo v = = ist, 


das Steigungsmaß oo hat; bei weiter fallendem v nimmt w ab, 
fir w—>0 geht v—>0. 
G, wird begrenzt durch K,, K,. Jede Lösungskurve in G, 


verläuft im Falle «=>, bei (w,»)=(0,0) beginnend, mit ab- 
nehmendem w zwischen den Kurven v= 0 und v steigend, 
erreicht da, wo v = = ist, das Steigungsmaß oo; bei weiter zu- 
nehmendem v steigt an; v erreicht sodann bei w = 0 in Maximum, 
bei weiter zunehmendem w fällt v; w = w(v) erreicht ein Maximum da, 
Wo v= = w nimmt sodann ab bei weiter fallendem v, bis 


fir v—»>0 auch w—> 0 geht. 
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Für « = 1 ist der. Verlauf in allen 4 Gebieten G, bisG, spiegel- 


- bildlich zur Geraden v = 0. 


Es bleibt noch das Gebiet VIII zu besprechen. Jede in VIII 
verlaufende Lösungskurve bleibt in jedem Teilgebiet w=W(W >1, 
beliebig) bezüglich der v-Werte beschränkt. Fürw—> 1:0—> — oo; 
für w—> © :0—>0. Keine Lösungskurve erreicht für endliches w 
den Wert v= 0. Alle Kurven in VIII verlaufen, von (w, v) = (1, — oo) 
kommend, steigend und streben für w—> oo dem Wert v = 0 zu. 

Zu betrachten sind noch die Grenzen zwischen (I, II, III), (IV), 
(V, VI, VII) und (VIII), das sind die Geraden v= 0 und w=1: 
v = 0 und w = 1 sind Lösungen von (3), wie man an der Parameter- 


darstellung a(w + A)wov; = 


Fr = z+1.9— w?z erkennt. 

Es sind damit alle im Endlichen verlaufenden Lösungskurven 
bekannt, die durch Punkte (w, v) gehen, durch die nur eine Lösungs- 
kurve geht. Bei den Punkten (w, v)=(1,0), (0,0), (—8, ee) —_ 
der letzte ist der Punkt B, vgl. z.B. Abb. 3 — kann nun noch die 
Möglichkeit bestehen, daß auch w=1,v=0) baw. w=(0, v = 0) 
bzw. (w =—f,v0= Lösungen von (3) sind. Das ist 


auch wirklich der Fall, wie durch Einsetzen in die Parameter- 
darstellung zu erkennen ist. 
Im Unendlichen der (w,v)-Ebene geiegene Lösungen. Durch die 


Substitution w = =, v= ; ergibt sich aus (3): 


dy ay” (1 + Bx) 
(4) 

dx - a? (1-2) 
Es ist zu sehen, daß da unstetig ist, wo r= 0, r= 1, 
y = (1 — zp ist da unstetig, wo y = 0, 
Folglich können höchstens durch die Punkte (z,y) = (0,0), (1,0), 


(+ in der z,y-Ebene mehrere Integral- 


kurven hindurchgehen. Festzustellen sind hier nur diejenigen 
Lösungen y='y(z), die in ihrem Erklärungsbereich und Werte- 
vorrat z= 0 bzw. y = 0 enthalten, und die mit ihrer Lösungskurve 
die Geraden y= 0 und z= 0 nicht verlassen, da sie sonst schon 
bei der Betrachtung der Lösungen erfaßt sind, die im Endlichen 
der (w,v)-Ebene verlaufen. An (4) ist zu erkennen, daß die 
Geraden y= 0 und z = 0 Lösungen sind. Weitere Lösungen, welche 
die Geraden y = 0 und z = 0 nicht verlassen, können dann höchstens 
5* 
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noch die Punkte (c= 0, y = 0) und = 1, y= 0) sein, die auch 
beide die Parameterdarstellung von (4) erfüllen. 

“ "Im Unendlichen der (w,v)-Ebene und gleichzeitig im Unendlichen 
der (z, y-Ebene liegen noch die Punkte (w, y) = (0,0) und (z,r) = (0,0). 
Durch die Substitution » -; ergibt sich aus (8): 

dy _ ay (w+ pw 
zy—2*t 
daran ist zu sehen, daB (w = 0, y = 0) keine Lösung ist. Durch die 
Substitution w = = geht (3) über in: 
dv a(1 + 
dx — att! 
demnach ist ‘der Punkt (= 0, » = 0) eine Lösung. 

"Es ‘gibt also insgesamt folgende Integralpunkte und Typen von 
Integralluurgen der Differentialgleichung (3) [vgl. Abb. 3, in der die 
Verhältnisse für Zylinderwellen (d.h. k = 2, a = 1) dargestellt sind]: 
(w, )= = (0,0); ph die Schar der Int.-Kurven innerhalb G,; 
(w, 0) = (1 0); K. :" » » ” ” ” G, 3 


. {die Schar der Int.- : 
= (-B- { Kurven innerhalb } 
(z, y) = (0,0); K,; die Schar der Int.-Kurven innerhalb G;; 
(5) { = (1,0); ” » =6, ; 
(z, 0) = (0,0); (V)=G,; 
v=0; ” n ” (VII) =6, ; 

w=l1 

y=0; 

z=(. 


Abb. 3. Die Einteilung (5) in Gebiete G, für Zylinderwellen in Gasen 


IL Flüssigkeiten 
Für den Fall der Flüssigkeiten (n = — 1; 8 =— 2 bei Zylinder- 
wellen, #=— 8 bei Kugelwellen) liest man in gleicher Weise wie 
unter I unmittelbar aus (3) die Stellen ab, wo in der (w,v)-Ebene das 
Steigungsmaß dv/dw positiv, negativ, Null oder oo ist. 
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Unbestimmt ist nach (3) das SteigungsmaB in den Punkten 
(w, 0) = (0,0); (1,0); (- 
den einzigen Punkten, bei denen zugleich und un- 
stetig sind, so daB durch diese Punkte und nur durch sie. mehrere 

Lösungen durchgehen können. 

In $ 3 der oben genannten Arbeit von Prof. Bechert ist die 
Lösung von (1) für den Fall n = — 1 hergeleitet" worden. Sie lautet, 
umgeschrieben auf die Größen v und w: 


(6). =0; 


C ist eine Integrationskonstante. 

’ Gestiitzt auf die entsprechenden Überlegungen bei der Kurven-- 
diskussion für den Fall der Gase und unter Verwendung von (6) 
läßt sich leicht folgendes beweisen [es ist k = 3 (Kugelwellen) an- 
genommen; das Kurvenbild für k = 2 (Zylinderwellen) gewinnt man 
durch Spiegelung der Ebene z < 0 an der Geraden v = 0]: 

»= 0 und z=0 sind Lösungen von 5 — In dem Punkt 


B = 
mündeu genau 4 Integralkurvenäste K,, K,, K,, K, K, stellt 
für 0<2<k—1 eine monoton wachsende Funktion v = v(w) dar; 
k* 
1" 
K, ist fir 0< 2<k—1 eine monoton-fallende Funktion v(w), die 
für w—+»1 über alle Grenzen wächst und ständig oberhalb der 
Kurve = verläuft. K, ist fiir z>k—1 eine monoton wachsende 


0, vw— k) —> 


Funktion, die ständig oberhalb der Kurve v = = verläuft und deren 


v-Wert für w—> oo über alle Grenzen wächst. K, stellt firz>k—1 
eine monoton fallende Funktion dar; v—»0 für w— > oo. — Im 
Punkt (w, v) = (0, 0) münden außer der Geraden v = 0 genau 2 Lösungs- 
kurvenäste K, und K,. K, ist fir 0<w<1 eine monoton wachsende 


Funktion, deren v-Wert für w—1 über alle Grenzen wächst. 


K, verläuft ständig oberhalb der Kurve v = =. K, stellt fir w< 0 
eine von (w,v) = (0,0) an monoton wachsende Funktion v(w) dar, die 
ständig oberhalb v = Fe verläuft und für w—> —oo über alle 
Grenzen wächst. 3 
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Die Geraden v= 0, 2=0 und die Kurven K, bis K, teilen 
die (w,v)-Ebene in die aus der folgenden Abb. 4 ersichtlichen Teil- 
gebiete G, bis G, ein. Die innerhalb eines der Teilgebiete G, bis G, 
verlaufenden Integralkurven verlassen das Teilgebiet nicht [abgesehen 
von der Fortsetzungsmöglichkeit über den Punkt (w,v) = (1,0) hinaus]. 


v 
% 
4 4 4 
6 | 4 


Abb. 4. Die Einteilung der (v,w)-Ebene in Gebiete G, für n =— 1 (Flüssig- 
keiten) und k = 3 (Kugelwellen), schematisch. Für k = 2 (Zylinderwellen) ist 
die Ebene z < 0 an der Geraden v = 0 zu spiegeln 


G,. Jede Integralkurve innerhalb G, ‘vermittelt eine eindeutige 
Beziehung v = v(w). Von (w,v) = (—00,0) kommend verläuft die 
Integralkurve fallend, hat ein Minimum bei w = 0, steigt an und 
nimmt fürw—>1 den Wert v= 0 an. —G,. Durch jede Integral- 
kurve innerhalb G, ist eine eindeutige Besiehung w= w(v) gegeben. 
Von (w,v) = (1,0) kommend verläuft die Kurve (v-Achse als Abszisse, 
z-Achse als Ordinate gesehen) steigend, hat da ein Maximum, 


wov= = ist, fällt, verläuft dabei ständig oberhalb der Kurve v= > 
[4 
und hat für v—> oo den Grenzwert z=0. — G,. Die Integral- 


kurven erklären eindeutige Funktionen v = v(w). v —» oo für z— 0 
und für z——»oo. Die Kurven verlaufen ständig — der 


Kurve v = =. v hat ein Minimum bei w=k; dort ist v> a = 

G,. Jede Integralkurve vermittelt eine eindeutige Beziehung w = w(v). 
w—> oo fir v— 0 und für v—> oo. w = w(v) hat ein Minimum 
da, wou = =. — G,. Die Integralkurven stellen eine eindeutige Be- 
ziehung v = v(w) dar. v—> 0 für z— 0 und für 2— > oo. v(w) hat 
ein Maximum bei w = k. —G,. Durch die Lösungskurven in G, ist 
eine eindeutige Beziehung v = v(w) hergestellt. v—>0 für 2—>0 
und für z—> oo. v(w) hat ein Minimum bei w = k. —G,. Durch 
die Lösungskurven innerhalb G, ist eine eindeutige Beziehung w= w(v) 
vermittelt. z—->0 für v—0 und fir v—>oo. w=w(o) hat 
ein Minimum da, wo v = =. — G,. Die Lésungskurven stellen ein- 
deutige Funktionen v = v(w) dar. v—+»oo für 2—>0 und für 
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z—> —oo. v(w) hat ein Minimum bei w = 0. —G,. Die Lösungs- 
kurven innerhalb G, geben eine eindeutige Beziehung w = w(v) 
w—> —oo fir v—>0 und fir w=—w(v) hat einen 
Höchstwert da, wo v = - 

Damit sind für den Fall n=— 1 sämtliche im Endlichen der 
(w,v)-Ebene verlaufenden Integralkurven von (3) zusammengestellt. 
Wie im Fall I geben im Endlichen der w, v-Ebene ferner die Integral- 


purikte (w= 0, v = 0); (w = 1, v = 0); ps) Vörungen. 


Im Unendlichen der (w,v)-Ebene liegen die gleichen Lösungen 
vor wie im- Fall I. 

Es ergibt sich somit das in Abb, 5 gegebene qualitative Bild 
vom Verlauf der Integralkurven von (8) in der (w, 2)-Ebene. 


| 


BSS 


Abb. 5. Der Verlauf der Integralkurven von (3) für n =— 1.(Flüssigkeiten), 
k = 3 (Kugelwellen). Bei k = 2: Spiegelung der Ebene z <0 . 
an der Geraden v = 0. Nur die ausgezogenen Linien sind Integralkurven 


Als Zusammenstellung aller Integralkurven und Integralpunkte 
können wir Gl. (5) auch für den Fall der Flüssigkeiten benutzen, 
wenn wir bedenken, daß jetzt 6 = — 'k ist, daß die Gebietseinteilung 
jetzt G, bis G, umfaßt und daß es 6 ausgezeichnete Kurven K, bis K, 
gibt. Die für die Flüssigkeiten geltende Bedeutung der Kurven K 
und der Gebiete G ist in Abb. 4 angegeben, sie ist anders als im 
Fall der Gase. | 

Es sei nun ein Überblick darüber gewonnen, welche Integrations- 
„konstante C in (6) den einzelnen Kurvenscharen der Abb. 5 zugehört. 
Hierzu verfolge ich den Wert von C = C(w,v) längs der zwei Geraden 


w=k; w= 0; und längs der Kurve of, durch die jede inner- 
halb der einzelnen Gebiete G, bis G, verlaufende Integralkurve 
insgesamt genau einmal geschnitten wird. Aus (6) ergibt sich: 
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Ferner gilt mit der Bezeichnung W = =; 


z 


k)+ 


dc a[d 
= v(k - 1)(v — kW} 
Ist die Differentiation längs der Kurve v = gemeint, sofolgt 
daraus: 


—2W fürk=2, 
(8) = 1(1 — — (1 — kt) = 12W? 


a) k = 2 (Zylinderwellen). Langs der Geraden w =k hat die 
Funktion C(w, v) = C (2, v) folgendes Verhalten gemäß (6) und (7) 
(2-25 =20-4): C(2; v) ist für 
(Punkt B!) eine monoton fallende, für v > 4 eine monoton steigende 
Funktion von v, die für v—> +00 den Wert oo, für v= 0 und 
a8 den Wert 0, fir v= 4 den Wert — 16 hat. — Langs 


der Geraden w=0 ist C(0,v) = vi. — Langs der Kurve — 


hat die Funktion C (w, )= C (w, W (w)) folgendes Verhalten gemäß (6) 
und (8): C (w, W (w)) ist in jedem einzelnen der Gebiete: —co < w< 0; 
0<w<1;1<w<2; 2<w< +00 eine monotone Funktion von w. 
An den Grenzen dieser Gebiete gilt C—» —oo bei w—» — oo; 
C=0 bei w= 0; — -o bei w=1; C= —16 bei w= 2; 
C—> 00 bei w—> +00. 

Daher stellen gemäß (6) die Kurven K,, K,, K,, K, die 


Funktion — 2%) = — 16 dar (es ist C = — 16), die Kurven K,, K, 
stellen die Funktion v = a dar (es ist C = 0). 

b) km 3 (Kugelwellen). Langs der w=3 zeigtC (w,v) = C(8,v) 
folgendes Verhalten gemäß (5) und (6) Sag = 3 = =) —3 =) 


- v) ist eine in den einzelnen Ge- 
bieten -o <v< <v< <v<o monotone, 
ah. 
C=0 bei v= ; C—> wfürv —> o.— 


hat C (0, v) den Wert v3, — Langs der Kurve v = = hat nach (6) und (8) 
C (w, v) = C (w, W (w)) folgendes Verhalten: C (w, W (w)) ist in jedem 


‘ 
+ 
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Teilgebit -—o <w<0;0 <w<1; 1<w<3; 3<w<o eine 
monotone Funktion von w. C —> oo für w— > +00 undfürw—> 1; 
C—+ 0 bei w—>0; —- bei w—> 3. 

Daher gehört zu dem Kurven K, bis K, die Integrations- 
konstante 3°/4 und gemäß (6) die Gleichung: o(e -3%) = =; 
die Kurven K,, K, stellen einen Zweig der Funktion v =3 Te = 0)dar. 


Diese Zuordnung zwisehen Integrationskonstanten und Integral- 
kurven sei nochmals in der folgenden Abb. 6 übersichtlich dargestellt. 


a) 


b) 


Abb. 6a.- Zuordnung der Konstantenwerte C (sie sind in der Hauptfigur an die 
Kurven geschrieben) zu den Integralkurven für n = — 1, k = 2 (Zylinderwellen 
in Flüssigkeiten). ‘Die kleinere Figur gibt die (v, w)-Werte der in der Hauptfigur 
hervorgehobenen Kurven 
Abb. 6b. Dasselbe wie in Abb. 6a, aber fürn =-—1,k =? 
(Kugelwellen in Flüssigkeiten) 


§ 2. Beschreibung der den Lösungen von (1) 
entsprechenden Wellenausbreitungen 
Die Gl. (lc), die bei Kenntnis des Zusammenhangs zwischen q 
und 9, also zwischen v und w, die Bestimmung der zugehörigen 


Werte von 7, &, ix? ermöglichen, lauten, auf Grund der Be- 


ziehungen [vgl. (1a)]: 


z 


| 
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umgeschrieben in v und w: 


1-a ? 


(9) n=Cexp.( aw); 


zv+Bwz 
(10) 
r -1 


Durch Einsetzen von (3) folgt aus (9): 


(12) n = C exp. (feau); 
Einsetzen von (12) in (10) gibt: 
(13) ger (frau); 


Einsetzen von (12) in (11) gibt: 


Es ist an (9), (10), (11) zu sehen, daß die Größen 7 und £, die aus 
einer bestimmten Lösung sich ergeben, für diese Lösung nie das 
Vorzeichen wechseln; +7 tut dies auch nicht, wenn für die be- 
treffende Lösung z immer < 0 oder immer > 0 ist, was für die in 
den einzelnen Gebieten G,, ... verlaufenden Lösungen zutrifft. § kann 
durch geeignete Wahl des Wurzelvorzeichens [vgl. (10)] stets als 
positiv angesehen werden. — Im folgenden wird dies nicht immer 


einzeln erwähnt, da es für alle zu besprechenden Lösungen gilt. 


I. Gase 


Zuerst. sei die Umdeutung der Lösungen von (3) auf die ent- 
sprechenden physikalischen Vorgänge für den Fall der Gase (n> 1; 
8 > 0) besprochen. Es werden nacheinander die in (5) aufgezählten 
Lösungstypen durchgegangen. Vgl. dazu auch Abb. 1, 2 und 3. 

Die 4 in (5) genannten Lösungsgeraden: v = 0; wel; y = 0; 
z= (0 geben keine physikalisch brauchbaren Lösungen, weil für 
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endliches 7 entweder oder ; überall Null oder oo wird. Von 
den Lösungspunkten ist aus dem iia Grund nur: 


wo) = (- 


brauchbar. Wegen (1a) erhält man in diesem Fall: y=-g=+1; 
das ist eine in $ 2 der Arbeit von Prof. Bechert behandelte Lösung. 

Es sind, wie dies auch in der Kurvendiskussion in § 1 
geschehen ist, nur solche Lösungen. der Differentialgleichung (3) zu 
betrachten, die eine Beziehung zwischen reellen v- und w-Werten 
darstellen. Denn es haben nur die Lösungen physikalische Be- 
deutung, die reelle Werte für 7? u, o liefern. Für n>1 
sind 7 und & bei reellem u, g reell, da Yn—1 reell ist. Wegen 

r -1 1 

sind daher dann nur Lösungen mit reellem z,v in Betracht zu 
ziehen. Ist n<1, so sind Yn—1,n,£ rein imaginär bei reellem u, o. 
Daher sind auch in diesem Fall wegen der im letzten Satz genannten 
Beziehungen nur Lösungen mit reellem v, z physikalisch zulässig. 

Für Gase liefern die Lösungen in den Gebieten @,,G, nach (10), 
(1b) imaginäre Dichten, wenn die Strömungsgeschwindigkeit reell ist. 
Diese Lösungen sind daher unbrauchbar. 

Die Lösungen im Gebiete G,. Zu. ihrer PORN auf den 


Zusammenhang zwischen den Größen n,$, ia betrachte man v 
als Parameter [weil w = w(v) eindeutig ist, und nicht v = v(w)], 
und suche zunächst 1, £, Zn in ihrer qualitativen Abhängigkeit. 
von v zu bestimmen, speziell fürv—> 0, v-——> oo; man untersuche 
das Vorzeichen der Funktionen und ihrer 1. Ableitungen nach ». 

n= C exp. ( f ® dw) I hat immer das Vorzeichen der Konstante C, 


da f @dw immer reell ist; die reelle oder komplexe Integrations- 
konstante steckt in C. Für alle v hat dn/dv umgekehrtes Vorzeichen 
zu C. 4 ist monotone Funktion von v; 4(v—» 0) endlich +0, 
n— co) > 0. & ist nach (10) für alle » entweder nur positiv 
oder nur je nach Wahl des Nach (13) 
hat E- =£ wir ~- gleiches Vorzeichen wie £, da 4 = — und der Nenner 


von W zugleich das Vorzeichen wechseln. & in eine monotone 
(überall oder überhaupt nicht 0 werdende) Funktion von v; § (v0 — 0) 
und §(v—-» oo) sind endlich. 
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i — ist nach (11) immer einerlei Vorzeichen. — Von v=0, 


Fr He oo herkommend, fällt der Absolutwert von iw bis zu 
einem Mindestwert, steigt an und strebt für v—>0o einem end- 
lichen von Null verschiedenen Grenzwert zu. 


Da 4 und £ monotone Funktionen von v sind, —— aber in 
einem gewissen Bereich fiir je zwei verschiedene o-Werte gleichen 
Wert hat, folgt, daß die Funktionen = und = 


t 

( zweideutig sind (vgl. Abb. 7); die 
ee Lésungen haben daher keine phy- 

7;  Sikalische Bedeutung, es sei denn, 

Bi a“ daB man den Bereich der Zwei- 

deutigkeit ausschlieBt und an das 

eindeutige Gebiet dieser Lösungen 


Abb. 7. Typ der Funktionen @ als Fortsetzung eine Lösung $ 7 
der Grundgleichungen anschlieBt, 


“ in einem gewissen Tey Bereich 


79 


in für die die eindeutig ist. Das ist grund- 
Integralschar des Gebiets G,; sätzlich möglich (vgl. z. B. die 
Gase, u< 0 entsprechenden Überlegungen zu 


den Lösungen im Gebiet G,). 

Die Lösungen im Gebiet G,. Die Überlegungen über die Lösungen 
im Gebiet G, lassen sich zum Teil unmittelbar auf die im Gebiet G, 
übertragen. n ist eine monotone Funktion von v; 7(0 —> 0) endlich, 
n(\v|—» 0) —> 0; dn/dv und haben gleiches Vorzeichen bei a = 1, 
‘verschiedenes Vorzeichen bei « = + & ist eine monotone, nie Null 
werdende Funktion von v; für o—»0O und |v|— oo bleibt £ 
endlich; d&/dv und & haben gleiches wenn 1, ver- 


schiedenes wenn « ist eine monotone Funktion 


von v; vnv=0 


oo 


i einem Grenzwert R zu, der = 0 ist, je nachdem, ob ae = +00 
oder —oo wird. 
Daher ist durch diese Lösungen für 
r 


ein eindeutiger Zusammenhang 
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erklärt; 7 und & sind monotone Funktionen von Tere und zwar 


steigen |7| und &(>0) mit wachsendem | | an; die beiden 


Funktionen sind bei 7 de nur für 12327 erklärt. 
—»0 für —> R| 


endlich, + 0, für ool 


7 


r 
& ist endlich, +0, und fur 


Diese nur für 2, R erklärten Funktionen 


nz ( 4 ) lassen eine Fortsetzung ihrer Erklärung bis —_ = 0 


t+t t+t 
zu, die den hydrodynamischen Grundgleichungen und der For- 
derung £ = §(n) genügt. Man setze dazu in dem Bereich zwischen 


| =0 und Te =|Rj (mit Einschluß der Grenzen): ' 
n=0; &=&(R), wobei §(R) den konstanten Wert bedeutet, den £ 


für ee = R in dem vorher betrachteten Gebiet Fesabe IR! an- 


nimmt. Man kann zeigen, daß bei —"— = R: 


t+t 
dn dg yee 
PER fir k = 3, 
dy Vn-1 dé wary: für k=2. 


Die Kurven = & (=): „= einer bestimmten Lösung - 
in G, haben daher, abgesehen 
von etwa noch auftretenden #7 


Wendepunkten, folgende Gestalt, ee 
wenn wir R>0, k = 2 wählen 


(vgl. Abb. 8). 4 ; 
Faßt man die horizontale R Fr 
Achse als r-Achse auf, so sind das 


zugleich die Kurven 1(r), &(r) bei Abb. 8. Typ der 
t+r=1; fir t+r+1 hat in Abhängigkeit von >-— für die 
man sich dann das ganze Kurven- _Integralschar in G,; Gase, R > 0, 

bild in Richtung der r-Achse im Zylinderwellen 
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Maßstab 1:(t + r) vergrößert zu denken. Über die Entstehung des 
darch eine solche Lösung beschriebenen physikalischen Vorgangs, 
d.h. für (¢+ t)< 0, gibt die Lösung keine Auskunft. Füri+r>0 
ist der Vorgang der: Zu einem bestimmten Zeitpunkt i+7=T 
liegt im Innern in 3) mit einem Radius RT 
um den Symmetrieort des physikalischen Vorgangs] keine Störung 
vor (7=0, &£=const). Nach außen hin nehmen Dichte und zentri- 
fugal gerichtete Strömungsgeschwindigkeit ‘monoton zu und streben 
für r —> oo einem endlichen Grenzwert zu. Durch das Nachaußen- 
streben der Teilchen nimmt die durchschnittliche Dichte in einem 
Zylinder bzw. einer Kugel mit einem Radius >RT ab; das ge- 
schieht in der charakteristischen Weise, daß das ungestörte Gebiet 
(in dem ja kleinstes £ und damit og ist) ständig auf Kosten von 
Teilchen des ‚gestörten Gebiets größer wird. Nach 82 der Arbeit 
von Prof. Bechert ist die Phasengeschwindigkeit eines bestimmten 


(n, &)-Wertes 
1 


Man kann zeigen, daß bei De gi En. somit läuft 
t ” 


Vk 
der Rand des ungestörten Gebiets von der Symmetrieachse bzw. dem 


weg, wenn $, die der Dichte g, im ungestörten Gebiet entsprechende 


Größe ist; ie — &, ist aber die zu 0, gehörige Schall- 


geschwindigkeit c, so daß also ee = ¢(@,). 

Durch R>O ist der ih genau in umgekehrter 
Richtung wie bei R>0 verlaufende Vorgang beschrieben. Er ist 
nur für {+r<0 erklärt; sobald die Wellenfront, d.h. die Grenze 
des ungestörten Gebiets, die Zylinderachse bzw. den Kugelmittel- 
punkt erreicht, beschreiben die unter der Voraussetzung § = £(n) 
abgeleiteten Gleichungen den weiteren Verlauf nicht mehr. 

| Die Lösungen im Gebiet G,. Die Integralkurven in G, vermitteln 
eine eindeutige Beziehung v = v(w). Man hat also jetzt 1, 
in Abhängigkeit von w zu untersuchen und findet: 

‘Far w—>0 strebt ir einem Grenzwert R zu; &(>0) und |n| 


sind mongton wachsende Funktionen. von. | 


Ter 


r 
die nur für 
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les 

<R<0, ‚wenn 7 <0, nur für > R> 0, wenn 7 >0, 

erklärt sind. Für —R bleibt n endlich, +0; § geht 

T gegen Null. Für at —> oo streben 7, £ gegen endliche Grenz- 
werte + ig 

ng 

setzung der Liésung: 0. 

= _ Der physikalische Vorgang ist für R>0 der: In einem Innen- 

wd gebiet ist keine Materie. Es schlieBt sich, mit der Dichte 0 beginnend, 

m ein gaserfülltes AuBengebiet an, dessen Dichte bei wachsender Ent- 

B- fernung von der Symmetrieachse bzw. dem Symmetriemittelpunkt 

et des Vorgangs zunimmt und für r—> co gegen einen endlichen 

oe Grenzwert strebt. Die Gasteilchen haben an der Grenze. von Innen- 

it und AuBengebiet eine nach außen gerichtete Gesehwindigkeit, die 

2. bei wachsendem r zunimmt und für r—> oo gegen einen endlichen 
Grenzwert strebt. Das gasfreie Gebiet vergrößert sich nun ständig, 
sein Rand entfernt sich von der Zylinderachse bzw. dem Kugel- 
mittelpunkt mit der Geschwindigkeit, die die Teilchen dort haben. 

ft Für R<0 liegt der raumzeitlich in umgekehrter Richtung ver- 
laufende Vorgang ¥ vor. 

m 

gegen zwei von Null "Werte R,. 

le Die Funktionen 1; -) und £ (4 pee ) sind mehrdeutig, -also 

ohne physikalische (wenn sie nicht geeignet eindeutig 
fortgesetzt werden). io sind in Abb. 9a dargestellt; es ist u > 0 

angenommen, daher 0. 

t Die Lösungen im Gebiet G,. Für v —>0, w —>0 geht 

e —R, für —o, w—>— oo geht 


Abb. 9a. Typ von q,u Abb. 9b. Typ von 9,u 
fir das Gebiet @,;' Gase für das Gebiet G,; Gase 
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R, und R, sind endlich nnd von Null verschieden. Die Lösungen 
geben mehrdeutige Funktionen 1 ( zz): für sie ergiht 


sich das in Abb. 9b dargestellte Verhalten, wobei wir offen lassen, 
ob die auf dem punktierten Stück auftretenden Stellen unendlich 
steilen Anstiegs von o, u vorliegen müssen. Bei einer Schar von 
Lösungen aus G, ist es der EM wie man zeigen kann. Für die 
Zeichnung ist u> 0, dahe 

Es bleiben noch die K za 

Lösung K,. K, verhält sich, von (w, v) = (0,0) kommend, wie 
eine Lösung in G,; es tritt wie bei den Lösungen in G, Mehr- 


deutigkeit der Funktionen 7 auf. Bei der An- 


näherung an den Punkt B=(w, v) = ( — f, Fr ) ändert sich 


das Verhalten; es gehen 2, 


r 
3 ry 
a gegen Null. Man bekommt den 
Verlauf der Abb. 10a. 
Lösungen K,,K,,K,. K, ver- 
halt sich, von (w, 0) = (— 00, 0) 
*“t kommend, wie eine Lösung in G, 
Abb. 10a. 9, u für die Integral- oder G,; K, von (w, |v]) = (— oo, 00) 
kurve K,; Gase herkommend wie eine Lösung 
in G,; K, von (w,v) = (0,0) her- 
kommend wie eine Lösung in G, oder G,. Beim Heranrücken an 
den Punkt B treten wieder erhebliche Abweichungen ein. Im 
ganzen gilt: 
Die Lösung K, liefert eine monotone Beziehung 


| 
die fir alle i =0 bzw. SO erklärt ist. Es ist &, 7 =0 
fürr=0, und &,7 endlich für r —> oo. 
Die Lösung K, liefert nur für | 7 2 >|R| erklärte Funk- 


tionen &, 7. Die Fortsetzung bis zu rd =0 erfolgt - ebenso wie 


bei den Lösungen in G,. Daher wie dort: Ungestörtes Innengebiet 
mit &,>0, n=0; sein Radius ändert sich mit der Geschwin- 
digkeit c(o,). Fir r— >» oo geht |) |—» 00. 

Die Lösung K, liefert nur für| 7 „|>1R | erklärte Funk- 


tionen &£,n. Fortsetzung wie bei den ‘en G,. Daher: Un- 
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gestörtes Innengebiet mit = 0; sein Radius ändert sich mit der 


1 
] = 0: 
Geschwindigkeit n(R)=u(R). Es gilt | n(R)|> 0, &(R)=0; 
für r—> co gehen &, |n|—> oo. Vgl. die Abb. 10b. 


Nebmen wir {+r>0, so haben wir in K, die Beschreibung 
einer Gasströmung, bei der im Symmetriemittelpunkt Dichte und 
Strömungsgeschwindigkeit Null sind; sie nehmen beide nach außen 
hin zu und streben mit wachsender Entfernung vom Symmetrie- 


Mm P? £ 
Bo) 7 


shy 


Abb. 10b. , 7 für die Integralkurven K, bis K,; Gase 


mittelpunkt gegen endliche Werte (die Strémungsgeschwindigkeit 
ist nach auBen gerichtet). Mit wachsendem ¢ sinken Dichte und 
Strömungsgeschwindigkeit an festem Ort; die Stellen fester Werte 
von Dichte und Strémungsgeschwindigkeit wandern mit zeitlich 
unveränderlicher, aber räumlich verschiedener Geschwindigkeit nach 


außen; jeder Punkt in Abb. 10b behält eben seinen Tey Wert bei. 


IL. Flüssigkeiten 


Wie die Diskussion des Verlaufs der Integralkurven, so ist auch 
die Umdeutung auf die physikalischen Vorgänge gegenüber dem Fall 
der Gase dadurch erleichtert, daß für n = — 1 die Lösung von (1) 
bekannt ‘ist: Gl. (6). Vgl. zum folgenden die Abb. 4—6. 

Wie schon erwähnt, sind zur physikalischen Deutung nur Lösungen 
brauchbar, die Beziehungen zwischen reellen w- und v-Werten dar- 
stellen. Da 7 und £ rein imaginär werden, wenn u und o reell sein 
sol'an, sind wegen (10) die Lösungen in den Gebieten G, und G, ur- 
brauchbar, denn sie führen bei imaginären, 7 auf reelles £ Es 
werden nun die in § 1, II aufgeführten Lösungen durchgesprochen. 
6 
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Von den Integralgeraden und Integralpunkten ist wieder nur die Punkt- 
lösung (w,v) = ß; brauchbar. Sie führt auf g=q=+1, 
und ist in § 2 der Arbeit von Prof. Bechert behandelt. 

Die Lösungen im Gebiet G,: | »| und | £| sind monoton wachsende 


r 
< R<0, 


r 
> R>O erklärt sind.. Bei 


gilt 7 —>0, £—> endlichen Wert; R ist der Wert von 


Funktionen von |. die bei > 0 nur für 


bei « < 0 nur für 


t+ı 
für w—» oc. Bei | — |—> gilt: 1, —» endliche Werte. 

Es ist zu beachten, daß nach (1b) g umgekehrt proportional 
zu & ist. Wie in $ 2,I können die Funktionen &, 4 über R hinaus 
fortgesetzt werden, hier durch: £ = const = §(R), 7 = 0. Es ist somit 
folgender Vorgang beschrieben (R > 0): An ein ungestörtes Innen- 
gebiet (9 = const; 7 = 0) schließt sich ein Außengebiet an, in dem 
nach außen hin die nach innen gerichtete Strémungsgeschwindigkeit 
zu- und die Dichte abnimmt. Dichte und Strémungsgeschwindigkeit 
streben für r — >» co endlichen Grenzwerten zu. Durch das Nach- 
Innen-Wandern der Teilchen wird das ungestörte Innengebiet (m 
dem größte Dichte herrscht) auf Kosten von Teilchen des Außen- 
gebiets ständig größer; sein Rand wandert mit der Geschwindigkeit R 
nach außen. — Bei R <0 liegt der raumzeitlich in umgekehrter 
Richtung verlaufende Vorgang vor. 


Abb. 11. Typ von 9,% für das Gebiet @,; Flüssigkeiten 


Die Lösungen im Gebiet G,. Hier sind zwei Fälle zu unter- 
scheiden. 1.v< =: |n| und |£| sind monoton steigende Funk- 


tionen von |, die nur für |R,|< 4 | < |R,| erklärt sind; 
r 


t+rt 


R, und R, sind endlich. Für R, gilt: 0, end- 


lichen Wert, Fir. —> R, gilt: |uj—>oo. R, R, 
sind >0, wenn u< 0, und umgekehrt: R,, R, < 0, wenn u > 0. 


R, ist der Wert von 


fir 2—> oo, R, entspricht 0; 
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2.0>k =: |n| und |&] sind monoton fallende Funktionen 


von | 74 |, die nur für |< erklärt sind, R,, R 
gehören zu den gleichen z-Werten wie unter 1. Auch die Zuordnung 
der Vorzeichen von R,, R, zu den Vorzeichen von % ist die gleiche 


wie vorher. Und es gilt wieder: 7 —»0, £—> endlichen Wert 


für R,, und |§|, |n|— oo für R,. 


Der Fall » =k = gibt unendlich vieldeutige Funktionen £, 7 
und scheidet daher aus. 

Die Lösungen im Gebiet G, sind nur in einem Teil ihres Be- 
reiches physikalisch möglich, weil siefür —— — R, in den beiden 
Fällen 1. und 2. für |&] und |»| unendlich große Werte liefern. 

Im Fall vsk = ist eine den Grundgleichungen gentigende 


Fortsetzung der Funktionen 7 ( ; = =) und £ ( z | über R, hinaus 
in der gleichen Weise wie bei den Lösungen in G, möglich: 
n=0, §= im für 


Auch über R, hinaus ist eine Fortsetzung der Funktionen £ und 4 
mathematisch möglich, physikalisch aber ohne Bedeutung, weil |£| 
und || bei R, unendlich werden. 


Für o<k“ und «<0 ist folgender Vorgang durch’ eine 


Lösung in G, beschrieben: An ein Innengebiet (r < R,(t + r)), in 
dem konstante Dichte und keine Strömung vorliegt, schließt sich ein 
Zwischengebiet (R,(¢ +r)<r<R,(t+r)) an, in dem nach außen 
hin die Dichte ab- und die Strömungsgeschwindigkeit, nach innen 
gerichtet, zunimmt. Hierdurch wächst der Radius des ungestörten 
Innengebiets bei gleichbleibender Dichte mit der Geschwindigkeit R, 
auf Kosten von Teilchen des Zwischengebiets. Dichte und Strömungs- 
geschwindigkeit würden im Zwischengebiet an dessen äußerem Rand 
theoretisch: die Grenzwerte 0 7 oo erreichen; natürlich gilt die 


Zustandsgleichung Pp =P — = dann längst nicht mehr, so daß das 


Verhalten der Lösung bei R, nur mathematisches Interesse hat. [Der 
äußere Rand wandert mathematisch mit der Geschwindigkeit R, > R, 

nach-auBen; über ihn weg strömt aus dem sich anschließenden Außen- 
gebiet Materie mit höchster (unendlicher) Geschwindigkeit und kleinster 
6* 
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(verschwindender) Dichte in das Innengebiet ein), — Der Vorgang 
ist durch die Lösung in G, nur für t + r) > 0 erklärt. Beiv<k = 


und bei wu > 0 ist der raumzeitlich in umgekehrter Richtung ver- 
laufende Vorgang beschrieben, und zwar nur bis zum Zeitpunkt 
(¢+ 7) = 0. 


Für v> on und u < 0 ist folgender Vorgang beschrieben: In 


einem Außengebiet (r > R,(t+r)) ist konstante Dichte und keine 
"Strömung. Es schließt sich ein Zwischengebiet (R,(t+r)<r<R, (t+7)) 
an, in dem nach innen zu Dichte und nach innen gerichtete Strömungs- 
geschwindigkeit ab- bzw. zunehmen. Im Außengebiet an dessen Rand 
wandern fortlaufend Teilchen nach innen: Der äußere Radius des 
Zwischengebiets wächst mit der Geschwindigkeit R,. Im Zwischen- 
gebiet haben an dessen inneren Rand Dichte und Strömungsgeschwin- 
digkeit die Grenzwerte 0 bzw. oo. [Über diesen inneren Rand, dessen 
Radius mit der Geschwindigkeit R, zunimmt, wandert nach dem 
mathematischen Bild Materie mit verschwindender Dichte und unend- 
licher Strömungsgeschwindigkeit in das Innengebiet (r< R, (¢+1))}. — 
Die Lösungen beschreiben diesen Vorgang nur vom Zeitpunkt 
(+n)=0 an. 


Bei v> = und > ist der raumzeitlich in umgekehrter 
z 


Richtung verlaufende Vorgang beschrieben, und zwar nur bis zum 
Zeitpunkt t + r) = 0. 
Die Lösungen im Gebiet G,. || und | €| sind monoton wachsende 


= r r 
Funktionen von | ee |, die fir alle ee =0, je nachdem 


Abb. 12. Typ von 9, u für das Gebiet G,; Flüssigkeiten 


ob u=0, erklärt sind. », & gehen gegen endliche Werte für 2—> oo; 
&—» 0, n-—— endlichen Wert für z—->0. Durch die Lösungen 
in G, wird folgender Vorgang beschrieben: Für u < 0 (d.h. wenn 
die Integrationskonstante C so gewählt ist, daß iC < 0): Die Dichte, 
die für r -—> oo einen endlichen Grenzwert + 0 hat, geht für r — 0 
gegen unendlich. Die nach innen gerichtete Strömungsgeschwindig- 
keit, die für r—> oo und r—>0 endliche Grenzwerte hat, nimmt 
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nach innen hin ab. Es wird ständig mehr Materie im Innern zu- 
' sammengedrängt. Trotz der im Innern unendlich werdenden Dichte 
befindet sich in einem Innengebiet mit endlichem Radius (kugel- 
“ symmetrischer Vorgang) bzw. endlichem Radius und endlicher Achsen- 
länge (Zylindersymmetrie) nur endlich viel Materie. Wegen der 
für r—> 0 endlich bleibenden Strémungsgeschwindigkeit verschwindet 
dort ständig Materie, und zwar endliche Massen in endlichen Zeiten. — 
Für u > 0 läuft der Vorgang raumzeitlich umgekehrt ab: es entsteht 
innen bei r = 0 Materie, die nach außen strömt. 

Die Lösungen im Gebiet G,. Sie verhalten sich für v—> 0 wie 
eine Lösung in G,, für 2—> oo wie eine Lösung in G,. 

Der Vorgang ist bei u< 0 der: In einem Innengebiet herrscht 
im Zentrum unendliche Dichte und endliche Strömungsgeschwindig- 
keit. Im Zentrum verschwinden in endlichen Zeiten endliche Massen. 
Nach außen hin wächst die nach innen gerichtete Strömungsgeschwin- 
digkeit bis zum Rand des Innengebiets hin mathematisch über alle 
Grenzen, und die Dichte nimmt bis zu 0 ab; vgl. hierzu die Be- 
merkungen zu den Lösungen im Gebiet G,. (Über den Rand des 
Außengebiets hinweg, der mit der Geschwindigkeit R vom Zentrum 
sich entfernt, wandert nach dem mathematischen Bild Materie mit 
unendlicher Strömungsgeschwindigkeit und verschwindender Dichte 
von außen nach innen). — Bei u > 0 raumzeitlich umgekehrt ge- 
richteter Vorgang. 

Die Lösungen im Gebiet G,. |n| ist eine monoton wachsende, 
|€| eine monoton fallende Funktion von | rss |; die nur fir 
lar 3 | > |R| erklärt ist. R ist der Wert von | er: | für |p| —> oo. 


Abb. 13. Typ von 9, u für das Gebiet G,; Flüssigkeiten 
£(00), n(00) sind endlich, +0; &(R) ist endlich, +0, n(R)= 0. 


R>0, wenn u>0; R<0, wenn u < 0. — Eine mögliche Weiter- 
erklärung der Funktionen £ und u über R hinaus: 


Es ist somit folgender Vorgang beschrieben (R>0): An ein 
ungestörtes Innengebiet ( = const, u = 0) schließt sich ein Außen- 
gebiet an, in dem Dichte und nach außen gerichtete Strömungs- 
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geschwindigkeit zunehmen und für r—> oo endliche Grenzwerte 
haben. Durch das Nach-außen-Wandern der Teilchen vergrößert ° 
sich der Radius des ungestörten Innengebiets, in dem kleinste Dichte 
ist, mit der Geschwindigkeit R. — Bei R< 0: Raumzeitlich in um- 
gekehrter Richtung verlaufender Vorgang. 

Die Lösungen im Gebiet G,. |n| ist eine monoton fallende 


Funktion von | die nur für IBl<| |R,} erklärt 
für z— > 0 baw. oo. 


t+r 


ist; R,, R, sind die Werte von 


n(R,) = 00, n(R,)= 0. |§| hat for 


fallt. bei wachsendem |; hat ein Minimum, steigt an und hat 


--» R, den Wert oo, 


für =. - » R, einen endlichen Wert +0. Es sind R,, R, 20, 


je nachdem ob u=0. In der Nähe von R, ist die Lösung offenbar 
physikalisch unbrauchbar. 

Es liegt folgender Vorgang vor: Über eine mit der Geschwindig- 
keit R, nach außen wandernde „Stelle“ r, = R, (t+ 1) strömt mit 
der Dichte 0 und der Geschwindigkeit oo Materie von innen nach 
außen. In einem Zwischengebiet steigt die Dichte an, hat ein 
Maximum und fällt bis zu einem von Null verschiedenen Wert wieder 
ab, während die Strömungsgeschwindigkeit monoton zum Wert 0 
abfällt. An das Zwischengebiet anschließend als mögliche Fort- 
setzung: ungestörtes Außengebiet (o = const, u = 0). Der zwischen 

rn=R(+r und 
"befindliche Wellenberg wandert nach außen, seine Breite r, — 1, 
wächst proportional mit der Zeit ((+r). — Bei R< 0: Raum- 
zeitlich umgekehrter Vorgang. 

Die Lösungen im Gebiet G,. Es ist RZ 0, je nachdem, ob u= 0. 
. Der nur angenähert physikalisch realisierbare Vorgang ist bei R > 0 
der: In einem Innengebiet herrscht im Zentrum unendliche Dichte 
und endliche, von Null verschiedene Strömungsgeschwindigkeit. Die 
Dichte nimmt bis zum äußeren Rand des Innengebiets, der mit 
konstänter Geschwindigkeit nach außen wandert, bis zum Wert 0 
ab, und die Strömungsgeschwindigkeit bis zum Wert oo zu. Über 
den Rand des Innengebiets strömt ständig Materie nach außen; im 
Zentrum entquillt ständig Materie. Bei R< 0: Umgekehrter Vorgang. 

Die Lösungen K,, K,, K,, K,. Im Fall der Zylinderwellen 
handelt es sich um Zweige der in § 4 der Arbeit von Prof. Bechert 
besprochenen Lösungen. Es sind solche, zu: denen in (6) die Inte- 
grationskonstante‘C = — 16 gehört (vgl. $ 1,ID; an Lösungskurven, 
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rte zu denen die Integrationskonstante C = — 16 gehört, gibt es außerdem 
ort ' noch eine im Gebiet G, und eine im Gebiet G,. 
hte K,. (Für u <0): Dichte und Strömungsgeschwindigkeit sind 
M- monotone Funktionen von r. Die Dichte ist unendlich für r—-> 0, 
sie verschwindet für r—» oo. Die Strömungsgeschwindigkeit, nach 
de innen gerichtet, hat für r —> 0 einen endlichen Wert. + 0, fir r —> co 
ärt wird sie unendlich. Im Zentrum verschwinden in endlichen Zeiten 
endliche Massen. Daraus, daß sich das Kurvenbild von g(r) z. B. 
00. nach doppelter Zeit (¢ + z) im Maßstab 1:2 vergrößert, kann sofort 
geschlossen werden, daß die gesamte zwischen r = 0 und r=co be- 
ui findliche Masse sich im Laufe der Zeit durch Zuwachs aus dem Un- 
at endlichen her vermehrt. Für «>0: Raumzeitlich umgekehrter. 
"Vorgang. 
0, K,; u<0. in einem Innengebiet (r< R(t+ r)) nehmen 
ar Strömungsgeschwindigkeit (nach innen gerichtet) bzw. Dichte von 
innen nach außen hin von 0 bis oo zu, bzw. von oo bis 0 ab; über 
ig- den mit der Geschwindigkeit R nach außen wandernden Rand des 
nit Innengebiets strömt Materie mit unendlicher Geschwindigkeit und 
ch verschwindender Dichte von außen. nach innen. Auf das Zentrum 
in zu staut sich immer mehr Materie auf. — Bei u > 0: Raumzeitlich 
ler umgekehrter Vorgang. 
0 K,. An ein ungestörtes Innengebiet (u = 0, o = const) schließt 
te sich ein AuBengebiet an, in dem für r—> oo die Strömungs- 
en geschwindigkeit unendlich und die Dichte 0 wird. Der Rand des 
Innengebiets wandert mit konstanter Geschwindigkeit nach außen. 
(u <0) Bei u> 0: Umgekehrter Vorgang. 
"1 K, Dichte und Strömungsgeschwindigkeit sind monoton ab- 
nehmende Funktionen von Für r—»oo sind Dichte und 
0. Strömungsgeschwindigkeit endlich, für r—> 0 wird die Dichte un- 
0 endlich und die Strémungsgeschwindigkeit 0; (u < 0). Man bekommt 
te wieder den Typ der Abb. 12. Für «> 0: Umgekehrter Vorgang. 
ie Lösungen K,, K,. ‘Sie bilden nach § 1,II die Integralkurve 
it v=k. = für z< 0. Die Lösungen sind physikalisch unbrauchbar, 
0 
er weil sie unendlich vieldeutige Funktionen £, 4 von 5; geben. 
m 
i! Damit sind alle Lösungen der hydrodynamischen Gleich agen 
»t untersucht, die Zylinder- oder Kugelwellen entsprechen, bei denen 


die Dichte o nur von der Strömungsgeschwindigkeit u abhängt, uid 
bei denen die Zustandsgleichung p = p, + £e gilt. 
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Physikalisch brauchbar sind im Fall der Gase die Lösungsscharen 
in den Gebieten G,,G,, und die Lösung K,. Für diese Lösungen 
werden Dichte und Strömungsgeschwindigkeit nicht unendlich; bei 
den letzten zwei Lösungstypen wird die Dichte an einer Stelle Null, 
die im Fall G, wandert, im Fall X, deshalb fest bleibt, weil sie 
ständig mit r = 0 zusammenfällt. Die übrigen Lösungen geben ent- 


weder Unendlichwerden von o und u für - pans —> oo oder Mehr- 
deutigkeit von oe und u in Abhängigkeit von = _ ; die Lösungen 
der letzteren Art müßten also über ihr Eindeutigkeitsgebiet hinaus 
geeignet fortgesetzt werden (was grundsätzlich möglich ist), um 
physikalisch deutbare Lösungen zu geben. 

Im Fall der Flüssigkeiten sind physikalisch brauchbar die Lösungs- 
scharen in den Gebieten G, und G,, als schematische Darstellung von 
Strömungsvorgängen sind außerdem brauchbar die Lösungsscharen 
in G, und die Lösung K,; bei den beiden zuletzt genannten Typen 
wird o beliebig groß für r—> 0. Bei den übrigen Lösungen wird » 
an einer Stelle Null, und u geht dort gegen oo; sie können also 
höchstens in einem Teil ihres Veränderlichenbereiches zur Dar- 
stellung physikalischer Vorgänge dienen. 


Gießen, Institut für theoretische Physik. 


(Eingegangen 27. November 1941) 
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